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ABSTRACT 


Mean  cross-shore  wave  height  transformation  and  alongshore  currents  observed 
on  near-planar  and  barred  beaches  are  compared  with  predictions  based  on  the  nearshore 
numerical  model  DelftSD.  DelftSD  solves  the  two-dimensional,  depth-averaged, 
momentum  balance  (2-DH)  between  forcing  (by  breaking  waves  and  variations  in  mean 
surface  elevation),  changes  in  momentum  flux,  bottom  stress  and  lateral  mixing.  The 
observations  were  acquired  on  the  near-planar  California  beaches  at  Torrey  Pines  and 
Santa  Barbara  and  the  barred  beach  at  Duck,  N.C.,  and  include  a  wide  range  of  conditions 
with  maximum  mean  currents  of  1.5  m/s.  The  model  has  two  free  parameters,  a  depth 
dependent  breaking  term,  y ,  and  the  bed  roughness  length,  .  An  empirical  formula  to 
determine  y  a  priori  from  the  deep-water  wave  steepness  and  bed  slope  is  developed, 
showing  good  agreement  in  the  wave  height  transformation.  Including  rollers  in  the  wave 
forcing  results  in  improved  predictions  of  the  observed  alongshore  current  structure  by 
shifting  the  predicted  velocity  maxima  shoreward  and  increasing  the  velocity  in  the 
trough  of  the  bar  compared  with  model  predictions  without  rollers.  On  near-planar 
beaches  and  high-energy  events  on  barred  beaches,  a  one-dimensional  (alongshore 
uniform  bathymetry)  model  performs  as  well  as  2-DH.  On  barred  beaches  under 
moderate  conditions  when  alongshore  non-uniform  bathymetry  prevails,  the  2-DH  model 
performs  better  than  the  1-D  model,  particularly  in  the  bar-trough  region.  Wave  forcing 
balances  the  bottom  stress  with  a  second  balance  between  alongshore  variation  in  the 
mean  surface  elevation  (pressure  gradients)  and  the  inertia  of  the  alongshore  current. 
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I.  INTRODUCTION 


A.  MOTIVATION 

On  planar  beaches,  as  depth  decreases,  waves  shoal  then  break  and,  depending  on 
the  slope  of  the  beach,  create  a  narrow  to  broad  cross-shore  distribution  of  breaking 
waves  known  as  the  surf  zone.  On  barred  beaches,  waves  break  over  the  bar,  reform  in 
the  trough  and  break  again  on  the  shore.  Resulting  nearshore  currents  generated  by 
obliquely  incident  breaking  waves  within  the  surf  zone  can  exceed  1  m/5  during  times  of 
large  incident  angle  or  high  wave  conditions.  Thus,  the  modeling  of  nearshore  circulation 
assumes  a  special  significance  for  military  operations,  where  the  success  or  failure  of  an 
amphibious  operation  or  infiltration/exfiltration  of  Special  Forces  is  directly  related  to 
our  ability  to  forecast  the  conditions  within  this  dynamic  area. 

Nearshore  flow  modeling  often  assumes  alongshore  uniformity  in  the  bathymetry. 
This  simplifies  the  momentum  equations  to  1-D  in  the  cross-shore  resulting  in  a 
mathematical  balance  between  the  wave  forcing  and  bottom  friction  in  the  alongshore 
and  the  wave  forcing  and  cross-shore  pressure  gradient  in  the  cross-shore  direction 
(Bowen  (1969),  Longuet  Higgins  (1970)  and  Thornton  (1970)).  Verifying  this 
assumption  using  a  1-D  model  in  conjunction  with  field  data  proves  surprisingly  accurate 
for  near-planar  beaches  (Thornton  and  Guza,  1986).  However  on  barred  beaches,  the  1-D 
assumption  has  proved  less  successful  when  compared  with  measured  data,  especially 
under  mild  wave  conditions  (Church  and  Thornton,  1993  and  Smith  et  ah,  1993). 
Random  wave  1-D  models  that  assume  wave  energy  dissipation  occurs  locally  at  wave 
breaking  (Battjes  and  Janssen,  1978;  Thornton  and  Guza,  1983)  result  in  local  radiation 
stress  gradients  and  alongshore  currents  predicted  over  the  bar  and  at  the  shore.  However, 
the  alongshore  current  profiles  obtained  during  the  Delilah  field  experiment  at  the  U.S. 
Army  Corp  of  Engineers  Field  Research  Facility  (FRF)  Duck,  NC  show  that  most  of  the 
time  the  current  maximum  was  within  the  trough  and  not  at  locations  where  intense  wave 
breaking  occurred. 

Varying  hypotheses  have  been  put  forth  to  explain  why  the  current  maximum 
occurred  over  the  trough.  These  include  rollers  (turbulent  layer  riding  on  the  forward  face 
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of  the  wave  erest)  during  breaking,  alongshore  pressure  gradients  and  mixing  of 
momentum  by  various  meehanisms.  Svendson  (1984)  introdueed  the  eoncept  of  a  surface 
roller  that  allows  wave  energy  to  be  transferred  to  a  roller  prior  to  dissipating.  The  time 
required  for  the  roller  to  dissipate  the  energy  creates  a  spatial  lag  between  wave  breaking 
and  the  energy  dissipation  (Naim  et  ah,  1990).  Comparisons  of  1-D  roller  models  with 
both  lab  (Reniers  and  Battjes,  1997)  and  field  data  (Lippmann,  Brookins  and  Thornton, 
1996)  have  shown  the  effect  of  the  roller  is  to  shift  the  wave  energy  dissipation 
shoreward.  The  roller  dissipation  mechanism  has  resulted  in  an  improvement  in  the  cross¬ 
shore  profile  of  the  alongshore  current  on  barred  beaches,  but  does  not  completely 
account  for  the  current  maxima  residing  in  the  trough. 

Since  most  beaches  fail  to  maintain  alongshore  uniformity  in  their  bathymetry  for 
any  length  of  time,  assuming  alongshore  uniformity  can  prove  to  be  too  simplistic. 
Therefore,  in  general,  the  full  horizontal  two-dimensional  depth-averaged  (2-DH) 
equations  should  be  used  to  compute  the  spatial  distribution  of  the  wave  and  flow  field. 
Putrevu  et  ah,  (1995)  showed  analytically  that  small  alongshore  variations  in  bottom 
topography  can  induce  alongshore  pressure  gradients  caused  by  alongshore  gradients  in 
breaking  wave  height  and  concomitant  set-up  that  can  force  an  alongshore  current  in  the 
trough.  Sancho  et  ah,  (1995)  used  a  quasi-3D  model  with  no  roller  dissipation  over 
smoothed  bathymetry  simulating  a  barred  beach  (Duck94)  to  examine  contributions  of 
the  different  terms  in  the  momentum  equations  to  the  flow  field.  They  found  that  the 
pressure  gradient  plays  a  significant  role  in  the  alongshore  flow  pattern,  resulting  in 
single-peak  alongshore  current  profiles.  Slinn  et  ah,  (2000)  also  examined  contribution  to 
the  momentum  equations  using  a  2-DH  depth  averaged  model  with  no  roller  dissipation 
over  smoothed  sinuous  bathymetry  designed  to  approximate  a  barred  beach  from  the 
Delilah  experiment.  Again  a  correlation  was  found  between  the  pressure  gradient  and 
flow  field  structure.  Neither  compared  wave  and  current  model  output  with  measured 
data  from  the  beaches  they  were  simulating. 

Currently  the  U.S.  Navy  utilizes  the  Navy  Standard  Surfzone  Model  (NSSM), 

which  provides  predicted  wave  heights  and  mean  currents  across  the  surf  zone.  The 

present  surf  model  is  limited  by  simplified  hydrodynamics  (1-D)  and  the  assumption  of 

alongshore  uniform  bathymetry.  Depth-averaged  (2-DH)  nearshore  numerical  models, 
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which  allow  for  the  inclusion  of  measured  bathymetry  and  complete  physics  in  the  cross¬ 
shore  and  alongshore  momentum  equations,  have  been  in  development  in  Europe  over 
the  past  two  decades  under  the  Marine  and  Science  Technology  (MAST)  program.  A 
number  of  mature  models  have  evolved  including  those  by  Delft  Hydraulics  (DelftSD), 
Danish  Hydraulic  Institute,  Service  Technique  Central  des  Ports  Maritimes  et  des  Voies 
Navigables,  HR  Wallingford,  Civil  Engineering  Department  of  the  University  of 
Eiverpool  (Nicholson  et  ah,  1997)  and  the  University  of  Delaware.  Model  differences 
reside  in  the  content  of  the  wave  driver  and  hydrodynamics.  Wave  drivers  vary  by  linear 
or  non-linear  monochromatic  or  random  waves,  wave  diffraction  and  wave  friction. 
Horizontal  circulation  models  solve  the  depth-averaged  Navier-Stokes  equations  using 
eddy  viscosity  turbulence  closure,  but  differ  in  their  evaluation  of  eddy  viscosity  using 
either  a  constant  or  allowing  spatially  variability.  Bottom  friction  in  the  current  models  is 
characterized  either  as  a  parameterization  of  bed  shear  stress  utilizing  wave-current 
interaction  or  separately  by  wave  or  current  forcing.  The  Naval  Postgraduate  School, 
Department  of  Oceanography,  Nearshore  Group  was  funded  by  the  Office  of  Naval 
Research  to  assess  the  DelftSD  model  as  an  upgrade  to  the  Navy  Standard  Surf  Model. 


B,  OBJECTIVES 

The  objectives  of  this  dissertation  are  to  assess  and  improve  the  ability  to  forecast 
wave  and  hydrodynamic  conditions  in  the  nearshore  through  the  evaluation  and 
improvement  of  the  nearshore  model  DelftSD.  DelftSD  is  an  improvement  over  NSSM  in 
that  it  includes  the  use  of  non-uniform  bathymetry,  an  improved  wave  driver  and  the 
depth  averaged  2-DH  momentum  equations.  Assessment  of  the  robustness  and  accuracy 
of  the  model  is  made  by  comparing  it  with  existing  field  data  from  both  near-planar  and 
barred  beaches,  investigating  the  gains  realized  by  advancing  from  1-D  to  2-DH 
hydrodynamics  with  3-D  bathymetry  and  assessing  the  implementation  of  a  roller  wave 
dissipation  mechanism.  Once  the  prediction  skill  of  the  DelftSD  model  is  established,  the 
focus  is  to  improve  our  understanding  of  the  contributions  of  the  terms  in  the  cross-shore 
and  alongshore  momentum  equations  to  the  nearshore  flow  field. 
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The  work  is  organized  as  follows:  Model  equations,  sensitivity  of  the  model  to  the 
free  parameters  and  the  model  formulation  are  presented  in  Chapter  II.  Chapter  III 
discusses  beach  types,  instrument  layout  and  wave/current  conditions  during  the 
experiments  conducted  at  Torrey  Pines  and  Santa  Barbara,  California  and  Duck,  North 
Carolina.  The  1-D  versus  2-DH  model  experiments  of  the  contributions  of  the  momentum 
terms  to  the  flow  field  are  presented  in  Chapter  IV.  Conclusions  and  recommendations 
are  presented  in  Chapter  V. 
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II.  MODEL  EQUATIONS 


A.  DELFT3D  EQUATIONS  (NO  ROLLER) 
1,  Wave  Model 


Assuming  steady  state  conditions,  no  mean  current,  no  wind  forcing  and  initial 
offshore  waves  described  in  terms  of  a  Raleigh  wave  height  distribution,  a  2-D  version  of 
Battjes  and  Janssen  (1978)  wave  transformation  model  is  used  to  solve  the  wave  energy 
balance  as  described  by  HISWA  (Holthuijsen  et  ah,  1989)  with  local  energy  dissipation 
described  as  a  linear  bore 


A 

dx 


=  -D,. 


-B 


(1) 


where  is  a  function  of  wave  direction  and  parameterized  as  a  frequency  spectrum  of 


fixed  shape,  and  represent  rectilinear  propagation,  represents  refraction,  B 
is  a  coefficient  of  0(1),  /  the  average  frequency,  the  depth  limited  wave  height 

and  E^^i  =  is  proportional  to  the  wave  energy  integrated  over  all  directions.  The 


fraction  of  breaking  waves,  Q^,  is  assumed  to  occur  as  a  delta  function  at  =ffmax 
a  truncated  Raleigh  distribution  which  results  in  the  implicit  relation 

(l-g>)  ,  (2) 

In  (a)  ffL 


where  is  the  maximum  wave  height  given  by 


0.88 


tanh 


/kh 


(3) 


in  which  k  is  the  wave  number  corresponding  to  /  .  In  shallow  water,  the  coefficient  y 


represents  the  ratio  of  the  breaking  wave  height  to  water  depth,  h,  and  is  the  only 
unknown  coefficient  in  the  model.  The  input  to  the  model  is  the  measured  bathymetry 
and  measured  wave  conditions  offshore.  The  measured  wave  conditions  include  the 
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significant  wave  height,  peak  wave  period,  mean  incident  wave  angle  and  the  directional 
spread. 


2,  Flow  Model 

The  fiow  module  solves  the  depth  averaged,  non-steady,  x-  and  y-  momentum 
equations  given  by 

du.  du.  drj  1  2  /  =  1,2 

dt  dxj  dx-  hp\^  dx-  j  7  =k2 

where  i,j  refer  to  horizontal  components  x,  y  ,  w.the  depth-averaged  velocity,  rj 

represents  the  surface  elevation,  andr.  the  combined  wave  and  current  bottom  shear 

stress  in  the  cross-shore  and  alongshore  direction.  The  radiation  stresses  in  both  the  cross¬ 
shore  and  alongshore  are  not  solved  directly,  but  formulated  as  proportional  to  the 
dissipation  divided  by  the  wave  celerity  (Dingemans  et  ah,  1987) 

55'.,.  D 

— (5) 

dXj  CO 

whieh  is  exact  in  the  alongshore  but  approximate  in  the  cross-shore.  This  approximation 
prevents  spurious  currents  from  forming  in  the  fiow  field.  The  eddy  viscosity, 

associated  with  lateral  mixing  is  held  constant  at  0.1  js .  The  bottom  shear  stress  is 
computed  using  Soulsby  et  ah,  (1993)  linearized  summation  of  combined  wave  and 


current  stresses 


r,  =/c(r„+r„) 


where  /c  is  a  fitting  function  of  the  Fredsoe  (1984)  non-linear  stress  model.  The  bed 

I 

shear  stress  due  to  the  current  is  calculated  using  a  quadratic  formulation,  r,,  =  p  Cf  \U\  , 

where  U  is  the  mean  current.  The  bed  shear  stress  coefficient  is  deseribed  after  White- 
Colebrook 


ISlog.o  ^ 


where  the  Nikuradse  roughness  length,  ,  typieally  related  to  grain  size,  is  used  as  the 
fitting  parameter  for  the  eireulation  portion  of  the  model.  The  bed  stress  due  to  the  waves 

-  /V 

is  modeled  as  =0.5 with  being  the  wave  frietion  faetor  and  f/^^^the 
amplitude  of  the  near-bottom  wave  orbital  veloeity  in  the  direetion  of  wave  propagation 
based  on  the  . 


B,  DELFT3D  EQUATIONS  (ROLLER) 
1.  Wave  Model 


The  inclusion  of  the  roller  dissipation  mechanism  requires  modification  to  the 
wave  transformation  model  (Reniers,  1999).  The  balance  for  the  short  wave  energy,  E^, 
to  describe  the  propagation  of  the  short  wave  groups  on  variable  bathymetry  is  given  bv 

dE^C^  cos(^)  dE„C^  sin(^)  _ 
dy 

where  0  is  the  angle  of  incidence  with  respect  to  the  x-axis,  x  the  cross-shore  distance 
and  y  the  alongshore  distance.  and  0  are  obtained  from  the  pre-computed  wave 

refraction  over  the  bottom  bathymetry  utilizing  HISWA.  The  wave  energy  dissipation, 
,  serves  as  input  in  the  balance  for  the  roller  energy  balance 

(9) 


dlEc  cos{0)  dlEc  sm{0)  ^  ^ 

^  ^  =  -D  +  £) 


dx  dy 

where  E^  is  the  roller  energy,  cthe  phase  speed,  =ct,  represents  the  roller  energy 
dissipation  expressed  by  (Deigaard,  1993)  and  z,  represents  the  wave-averaged  shear 

stress  between  the  roller  and  the  wave  interface.  The  shear  stress  for  a  steady  roller  is 
given  by  (Duncan,  1981) 

pgAsmiP) 


T,  = 


L 


(10) 


where  p  is  the  water  density,  g  is  gravitational  acceleration,  fi  represents  the  slope  of 
the  wave  front,  A  is  the  roller  area  and  L  is  the  wave  length  corresponding  to  the  peak 
period.  Roller  energy  area  is  related  to  the  roller  energy  by  (Svendsen,  1984) 
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(11) 


combining  equations  (2)  and  (3)  using  the  roller  energy  dissipation  results  in 

_  2gsin(yg)£^  (12) 

r 

c 

with  P  an  additional  free  parameter,  held  eonstant  at  0.04  throughout  the  assessment, 
eonsistent  with  values  found  by  Walstra  et  ah,  1996.  In  general,  adveetion  length  of  the 
roller  inereases  with  deereasing  p  .  E^,\s  taken  to  be  zero  at  the  upwave  boundaries. 

2,  Flow  Model 


The  momentum  equations  with  the  roller  ineluded  differ  in  two  ways.  First,  the 
radiation  stress  terms  inelude  a  roller  eontribution  (Renters,  1999) 


c 


k.k  E 

r  1 

2^-1 

_  k^  2 

c 

i  =  l,2 
^  ’  7  =  1,2 


(13) 


where  is  the  wave  energy,  k^kjjk^  axQ  the  direetion  eosines  in  both  the  eross-shore 
and  alongshore,  k  is  the  wave  number,  c  is  the  phase  speed,  is  group  veloeity  and  the 
eddy  viseosity.  Seeondly,  eddy  viseosity  eoeffieient  ,  assoeiated  with  lateral  mixing  is 
determined  using  a  formulation  by  Battjes  (1975) 


Z'  r\  \ 


o,=h 


D 


P  ) 


(14) 


where  the  roller  energy  dissipation,  D^,  is  obtained  from  the  roller  transformation. 
Hereafter,  the  eombined  eddy  viseosity  and  eddy  viseosity  eoeffieient  are  termed 
turbulent  mixing  ( ). 


C.  MODEL  FORMULATION 

The  numerieal  model  is  deseribed  in  detail  by  Stelling  (1984)  and  Verboom  and 
Slob  (1984).  It  is  a  finite  differenee  model  utilizing  an  Alternating  Direetion  Implieit 
(ADI)  teehnique  on  a  staggered  grid.  The  ADI  seheme  implieitly  solves  the  water  levels 

and  veloeities  in  the  x  -direetion  in  the  first  half  of  the  time  step  and  the  y  direeted  terms 
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in  the  second.  The  grid  spacing  is  varied  based  on  the  bathymetry  feature  needing  to  be 
resolved.  In  general,  10  grid  points  are  needed  to  capture  the  feature  of  interest.  For 
Delilah  a  grid  spacing  of  Ax  =3  m,  Ay  =6  m  is  used  to  capture  the  non-uniform 
bathymetry  in  the  cross-shore  and  alongshore.  For  Santa  Barbara  and  Torrey  Pines  a  finer 
resolution  was  required  ( Ax=l  m,  Ay  =2)  due  to  the  steepness  of  the  beach  in  the  cross¬ 
shore  and  the  longer  period  incident  swells.  The  time  step  utilized  {At=  1.5  s)  is  based 
on  the  Courant  number,  CFL,  for  wave  propagation 

CFL  =  2  At 

where  Ax  and  Ay  are  the  smallest  grid  spaces  in  the  physical  space. 

Model  bathymetry  is  created  from  bathymetry  measurements  taken  during  the 
various  field  experiments  using  a  triangular  interpolation  scheme.  The  criteria  for 
offshore  boundary  location  consists  of  minimizing  the  number  of  grid  points  required  to 
fill  the  domain  of  interest  and  selecting  a  water  depth  well  outside  wave  breaking.  This 
resulted  in  an  offshore  boundary  water  depth  of  4  m.  for  all  days  examined  except  the 
large  wave  days  at  Delilah  (Oct.  11-13)  where  wave  shoaling  at  h  =  A  m  required  the 
boundaries  be  moved  to  8  m.  The  daily  bathymetry  collected  during  Delilah  extended  to 
approximately  4  m  depth  with  a  pre  and  post  experiment  extension  to  8  m.  In  order  to 
develop  model  bathymetry  to  8  m,  a  composite  of  the  daily  measured  bathymetry  was 
merged  with  the  post  experiment  8  m  set.  This  proved  adequate  as  the  bathymetry  from  4 
m  depth  to  8  m  changed  little  between  the  pre  and  post  8  m  bathymetry  set.  A  smoothing 
scheme  is  implemented  over  the  model  bathymetry  for  days  when  the  bar  was  attached  to 
the  shore  (Delilah:  Oct.  7-9)  to  prevent  large  alongshore  depth  gradients  from  causing 
unrealistic  flow  patterns. 

The  boundary  conditions  for  the  wave  model  are  separated  into  those  required  to 
satisfy  the  HISWA  wave  model  (Equation  1)  and  those  for  the  wave  model  using  a  roller 
dissipation  mechanism  (Equations  8  and  9).  The  HISWA  offshore  boundary  requires 
input  of  the  significant  wave  height,  incident  wave  angle,  peak  wave  period  and 
directional  spreading.  Values  for  the  offshore  boundary  wave  input  are  derived  from  the 
8m  arrays  and  linearly  shoaled  and  refracted  to  the  4  m  water  depth.  The  wave  energy  at 
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the  lateral  boundaries  was  chosen  to  be  zero.  In  the  case  of  the  wave  model  that  included 
a  roller  dissipation  mechanism,  the  roller  energy  is  prescribed  zero  at  the  offshore  and 
lateral  boundaries. 

In  addition  to  prescribing  tidal  elevation  at  the  boundaries,  the  flow  model 
requires  two  boundary  conditions  for  closure  at  the  offshore  and  lateral  boundaries.  First 
a  Riemann  invariant  is  used  to  allow  disturbances,  created  by  differences  in  water 
elevation  between  the  boundary  and  inner  domain,  traveling  perpendicular  to  the 
boundaries  to  pass  in  a  non-reflective  manner.  The  second  flow  boundary  condition  is 
such  that  at  inflow  the  advection  of  momentum  in  the  tangential  direction  is  neglected.  At 
the  land-water  interface  it  is  assumed  no  flow  through  the  boundary  and  zero  tangential 
shear  stress  (free  slip). 

Steep  incident  wave  angles  0(20°)  can  result  in  large  setup  gradients  near  the 
lateral  boundaries  close  to  shore  resulting  in  unrealistic  circulations  and  flow  forcing  that 
can  propagate  into  the  model  domain.  To  prevent  boundary  generated  disturbances  from 
contaminating  the  flow  field,  the  model  domain  is  extended  in  the  alongshore  0(500)  m 
on  either  side  of  the  domain  area  of  interest.  The  extension  is  created  with  a  systematic 
increase  in  alongshore  grid  spacing,  culminating  in  a  Ay  of  50  m  at  the  lateral 
boundaries.  These  large  grid  steps  act  to  diffuse  the  water  level  gradients  and  prevent 
spurious  currents  from  propagating  into  the  domain.  The  bathymetry  in  this  portion  of  the 
grid  is  expanded  uniformly  alongshore  from  the  last  measured  cross-shore  profile  prior  to 
the  artificial  extension. 
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III.  FIELD  DATA 


A.  PLANAR  BEACH 

1.  Torrey  pines 

The  Torrey  Pines  Beach,  located  north  of  San  Diego,  California,  is  near-planar 
with  a  gentle  (1:50)  slope.  Due  to  shadowing  and  refraction  by  offshore  islands,  the 
incident  swell  direction  was  less  than  15  deg.  One -hour  records  were  analyzed  for  data 
acquired  on  Nov.  10,14,18,20-21,  1978  (17  hours).  The  peak  frequency  of  the  incident 
wave  spectra  varied  little  during  the  experiment  and  was  about  0.06-0.08  Hz  (7),=  12.5- 

16  s).  The  rms  wave  height  ranged  between  0.47  -  0.69  m  in  8  m  depth  and  initial 
offshore  rms  wave  height  (in  the  model)  was  set  equal  to  the  measured  hourly  value. 
During  this  experiment  and  all  other  experiments  described  herein  the  cross-shore 
transect  of  wave  heights  and  currents  were  measured  using  a  cross-shore  array  of 
pressure  sensors  and  Marsh-McBimey  electromagnetic  current  meters.  The  pressure  array 
extended  from  8  meters  depth  to  inside  the  surf  zone  with  the  current  meter  array  from 
2.5  meters  depth  to  inside  the  surf  zone  (Thornton  and  Guza,  1983).  The  data  were 
sampled  at  2  Hz.  The  tidal  range  is  2  m.  Only  data  over  the  high  tide  are  used  when 
instruments  are  submerged.  Surface  elevations  are  inferred  from  pressure  data  for  all  data 
sets  described  herein  by  applying  the  appropriate  linear  wave  theory  spectral 
transformation  (Guza  and  Thornton,  1980).  The  data  are  band-pass  filtered  from  0.05  to 

0.3  Hz  to  remove  high  frequency  noise  and  low  frequency  (surfbeat)  waves.  77^^,  =  VScr^ 

is  calculated,  where  cr^  is  the  surface  elevation  variance.  Bathymetry  profiles  were 

collected  daily  by  rod  and  level  survey  along  the  measurement  transect.  Table  4.1 
summarizes  wave  conditions  and  bottom  slopes  for  the  beaches  analyzed. 

2.  Santa  Barbara 

Leadbetter  Beach  located  in  Santa  Barbara,  California  is  a  mild  (1:25)  sloping 
near-planar  beach  with  a  tidal  range  of  2  m.  Protection  by  the  offshore  Channel  Islands 
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and  Point  Conception  to  the  north  results  in  a  narrow  aperture  for  ineoming  oeean  swell 
(less  than  9  deg).  Pressure  and  velocity  data  were  acquired  at  14  loeations  in  a  cross¬ 
shore  array  starting  in  3  meters  (Thornton  and  Guza,  1986).  Data  were  sampled  at  2  Hz. 
The  measurements  were  taken  over  high  tide  when  all  instruments  were  submerged.  The 
3  days  of  data  (Feb.  4-6;  13  hours)  utilized  for  model  comparison  consist  of  narrow- 
banded  incident  waves  with  periods  ranging  between  12.0-14.3  s  and  relatively  small 
wave  heights  {H^^  =  0.29-0.65  m).  Daily  nearshore  bathymetry  profiles  were  measured 

by  rod  and  level  surveys  at  5  cross-shore  range  lines  covering  an  alongshore  extent  of 
200  meters. 

B,  BARRED  BEACH 

1.  Delilah 

The  beach  at  Duck,  North  Carolina  is  typically  a  barred  beach  with  a  mean  tidal 
range  of  approximately  1  m.  The  mean  foreshore  slope  of  the  beach  is  approximately 
0.08  (1:12)  and  1:35  offshore  from  the  bar.  Wave  heights  during  the  Delilah  experiment 
were  measured  using  a  eross-shore  array  of  9  pressure  sensors.  The  data  were  sampled  at 
8  Hz  continuously  from  Oct.  7-19  (300  hours).  The  area  of  the  beaeh  where  the 
instruments  were  deployed,  known  as  the  mini-grid  region,  was  surveyed  daily  along  20 
eross-shore  profile  lines  spaeed  25  m  apart  near  the  instruments  and  50  m  apart  elsewhere 
resulting  in  an  area  approximately  550  m  in  the  alongshore  and  300  m  in  the  eross-shore 
to  a  depth  of  4  m.  Offshore  directional  wave  spectra  were  measured  using  a  linear 
alongshore  array  of  10  pressure  sensors  in  8m  depth  (Long,  1994). 

The  waves  and  resulting  alongshore  eurrents  varied  eonsiderably  during  the 
experiment,  Oet.1-20.  Waves  and  currents  were  mild  from  Oet.  7-8.  On  Oet.  9  a  strong 
frontal  system  from  the  south  generated  broadband  waves  at  relatively  large  angles  ('-'40 
degrees  from  the  south  in  8  m  depth)  resulting  in  a  strong  alongshore  eurrent  to  the  north 
(1.5  m/s).  The  wave  event  began  to  subside  on  the  afternoon  of  the  11**^,  but  on  the 
evening  of  the  12**'  large  swell  waves  began  to  arrive  from  the  south  due  to  distant 
Hurricane  Lilli;  these  waves  grew  to  a  maximum  on  the  morning  of  the  13*  (Hs-'2.5m) 
then  gradually  diminished  to  a  wave  height  of  --I m  on  the  morning  of  the  14*.  Although 
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the  wave  heights  diminished,  the  alongshore  eurrents  remained  strong  (>  Im/s)  for  the 
remainder  of  the  experiment.  Wave  direetions  changed  considerably  the  last  4  days  with 
the  incident  wave  direction  swinging  from  south  to  north  and  then  south  again  on  the  16* 
causing  the  alongshore  currents  to  change  direction.  The  swell  became  bi-modal  on  the 
17*,  persisting  through  the  18*  with  the  alongshore  current  flowing  to  the  north.  On  the 
19*,  a  small  frontal  system  passed  to  the  north  generating  waves  from  the  north  and 
currents  to  the  south. 

2,  Duck94 

The  Duck94  nearshore  experiment  was  conducted  during  Sept,  and  Oct.  1994  at 
the  Field  Research  Facility;  Duck,  North  Carolina.  The  beach  during  the  experiment 
could  be  characterized  as  a  barred  beach  (30-120  m  offshore).  Directional  wave  spectra 
were  again  acquired  using  a  linear  array  of  10  pressure  sensors  in  8  m  depth. 
Additionally,  a  13  element  cross-shore  array  of  pressure  sensors  was  used  to  measure 
wave  heights  spanning  the  width  of  the  surf  zone  to  8  m  depth.  The  data  were  sampled  at 
2  Hz  (Elgar  et  ah,  1997) 

The  data  selected  for  analysis  are  from  a  24-hour  data  run  acquired  on  September 
16*.  The  16*  was  a  mild  day  =  0.28-0.33  m)  with  a  long  period  swell  {T^  =13.7  s). 

This  combination  of  conditions  resulted  in  wave  steepness  values  less  than  0.002.  Wave 
data  with  steepness  values  less  than  0.002  was  required  in  the  development  of  an 
empirical  formula  for  determining  the  wave  breaking  free  parameter,  y ,  a  priori. 
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IV.  MODEL  RESULTS 


A.  MODEL-DATA  COMPARISON  TECHNIQUES 


Comparison  is  made  between  the  measured  and  eomputed  results  of  both  the 
eross-shore  transformation  of  the  root  mean  square  ( rms )  wave  height,  ,  and  the 

mean  eross-shore  (14)  and  alongshore  (v^)  variation  of  the  alongshore  eurrent  to 


determine  model  performanee.  Model-data  eomparisons  are  evaluated  by  first  ealeulating 
the  relative  rms  errors  defined  by 


s 


rms 


W- 


rms{Model)  ^  rms{Measured) 


)■) 


(16) 


where  (  )  denotes  average  over  all  values,  P  is  model  predietion  and  O  is  the  observed 
quantity.  The  s  value  is  then  subtraeted  from  one,  generating  a  model  skill, 
skill  =  1  -  ^  .  Skill  value  represents  an  rms  differenee  between  model  predietions  and 
observations  normalized  by  a  known  result,  with  a  value  of  0  indieating  the  model  has  no 
skill,  and  the  value  1  indieating  there  is  no  differenee  between  model  and  measured  data. 
The  statistie  is  not  well  behaved  for  low  foreing  eonditions,  whieh  have  maximum 
observed  it  and  along  an  instrument  transeet  <  0.2  m/s.  Observations  for  that 

transeet  are  disregarded  during  that  hour.  The  model  skill  aids  in  quantifying  model 
performanee  and  is  used  to  determine  optimal  7  and  during  model  ealibration.  In 
addition,  when  eomparing  large  amounts  of  observed  and  eomputed  data,  a  linear 
eorrelation,  r  and  absolute  rms  differenee,  s ,  are  used 


£  = 


( p  -  o 

y  rms(Model)  ^  rms{Measured) 


f) 


(17) 


Absolute  differenee  is  not  affeeted  by  small  observed  data  values,  and  provides  a  mean 
quantitative  value  for  the  disparity  between  the  observations  and  eomputed  results. 
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B,  MODEL  PARAMETERS 

1.  Wave  Breaking  Parameter  ) 

a.  Introduction 

The  objectives  of  this  section  are  to  determine  model  sensitivity  to  the 
depth  dependent  wave  breaking  free  parameter  y,  investigate  an  existing  empirical 
formulation  for  the  determination  of  y  a  priori,  and  to  develop  a  new  formulation  that 
expands  the  predictive  range  of  y  to  more  beach  types. 

Several  surf  zone  wave  transformation  models  take  into  account  random 
waves  and  describe  breaking  as  a  linear  bore  with  a  local  probability  of  wave  breaking 
function  (Battjes  and  Janssen,  1978,  Thornton  and  Guza,  1983).  Verification  of  the 
Battjes  and  Janssen  random  wave-breaking  model  was  conducted  by  Battjes  and  Stive 
(1985)  (hereafter  denoted  BS85)  using  both  laboratory  and  field  data.  They  determined 
an  empirical  formula  for  the  single  wave  breaking  parameter  in  their  transformation 
model 

7^535  =0.5 +  0.4  tanh(335o)  (18) 

where  the  deepwater  wave  steepness,  Sq  =  H ,  with  being  the  deepwater 
root  mean  square  wave  height  and  =  g jin  ,  the  wavelength  at  the  peak 
frequency,  .  The  data  used  by  BS85  (see  Figure  4.1),  produced  no  values  for  less 

than  0.01,  a  region  typically  associated  with  longer  period  swell. 

b.  Observations 

Wave  data  acquired  on  the  near-planar  California  beaches  at  Torrey  Pines 
(1:50  slope)  and  Santa  Barbara  (1:25  slope)  as  part  of  the  Nearshore  Sediment  Transport 
Study  (NSTS)  experiments  and  the  barred  beach  at  Duck,  NC  during  the  Delilah  (Oct. 
1990)  and  Duck94  (Sept.-Oct.  1994)  field  experiments  conducted  at  the  U.S.  Army  Corps 
of  Engineers  Field  Research  Facility  (FRF)  are  used  to  expand  the  work  of  Battjes  and 
Stive  (1985).  A  2-DH  version  of  the  wave  model  and  longer  period  swell  with  wave 
steepness  values  ranging  from  0.0009-0.02  are  used  to  extend  the  BS85  calibration 
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formula  for  y .  These  data,  dominated  by  long  period  swell,  resulted  in  an  extensive 
range  of  Sq  values  less  than  0.01  (a  region  not  yet  evaluated  for  empirieal  dependenee  of 
y  on  )  to  inelude  data  with  values  less  than  0.002  (a  region  where  y  beeomes 
independent  of  s^).  Figure  4.1. 

c.  Model-Data  Comparison 

The  wave  model  was  run  to  optimize  the  wave  breaking  parameter,  y^^^ , 
based  on  iteratively  minimizing  the  relative  error,  (Equation  16)  between 

measurements  and  the  model  predieted  along  the  eross-shore  transeet,  under  the 
eonstraint  B=\.  The  best  fit  values  for  the  wave  breaking  parameter, ,  are  determined 
using  one  hour  means  from  the  354  hours  of  observed  H data  aequired  during  the 
Torrey  Pines  (17  hours),  Santa  Barbara  (13  hours),  Delilah  (300  hours)  and  Duok94  (24 
hours)  experiments.  Waves  measured  offshore  during  the  one-hour  observation  are  used 
as  input  to  the  wave  model.  The  resultant  H model  output  values  are  eompared  with 
measured  data  from  the  various  beaeh  types  with  differenees  between  the  eomputed  and 
measured  eompared  statistieally. 


The  values  are  determined  using  the  empirieal  formula.  Equation  1 8 
using  the  same  350  hours  of  data.  Values  for  are  derived  from  the  deep  water  wave 
height  determined  by  reverse  shoaling  and  refraeting  measured  from  the  farthest 

off-shore  pressure  sensor  (indieated  with  the  subseript  r ,  Equation  19)  in  the  eross-shore 
transeet,  assuming  a  plane  sloping  bed  offshore  and  using  linear  wave  theory 
transformation  at  the  peak  period. 


H 


rmsQ 


=  H„ 


(c  ^ 

gf- 

2 

^  eos  6^  ^ 

fcos^o ) 

(19) 


where  subseript  0  indieates  deepwater,  represents  the  group  speed  and  6  is  the  mean 
ineident  wave  angle. 
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As  typical  examples,  model  predietions  are  eompared  as  a  funetion 

of  eross-shore  distanee  with  measured  values  (*)  from  Torrey  Pines,  Santa  Barbara, 
Delilah  and  Duek94  (Figure  4.2).  The  bottom  profile  is  given  as  a  referenee.  Model 
results  are  shown  for  both  (solid  line)  and  (dashed  line)  with  skill  value 

ineluded  (legend  box.  Figure  4.2).  Both  the  Torrey  Pines  and  Santa  Barbara  observations 
show  a  broad  eross-shore  distribution  of  wave  breaking.  The  Delilah  measurements  show 
two  breaker  regions;  the  first,  just  seaward  of  the  bar,  is  somewhat  broad;  the  seeond,  at 
the  beaeh  faee,  is  somewhat  abrupt.  Smaller  ineident  wave  heights  in  the  Duek94  ease 
result  in  a  single  breaker  region  at  the  beaeh  faee. 

Good  agreement  is  found  between  measured  and  modeled  for  all 
data  using  7^^,  for  both  the  near-planar  NSTS  data  and  the  barred  profile  Delilah  and 
Duek94  data  with  skill  >0.95  .  In  the  Torrey  Pines  ease,  the  range  of  7^^,  over  the  17- 
hour  span  analyzed  is  0.25-0.40.  Correlation  plots  between  modeled  and  observed 
for  all  sensors  using  7^^^  result  in  good  fit  with  r  =  0.94  and  s  =  0.91  at  all  sensors  (top 

left  panel.  Figure  4.3).  For  Santa  Barbara,  r  =  0.98  and  5  =  0.93  (top  right  panel.  Figure 
4.3).  The  data  from  Duek  is  a  eombination  of  mean  hourly  values  of  spanning  300 
hours  during  Delilah  and  one  24-hour  period  during  Duok94  (ineluded  to  provide 
additional  values  less  than  0.002).  Comparing  measured  and  modeled  derived 
from  y^p,  resulted  in  r  =  0.98  with  5  =  0.93 ,  (bottom  panels.  Figure  4.3).  Correlation 
eoeffieient  for  the  24-hour  segment  of  Duok94  data  is  0.72  for  the  y^^,  ease,  but  has  an 

equivalent  skill  to  Delilah.  Lower  eorrelation  eoeffieient  values  for  Duok94  result  from 
little  spread  in  the  wave  height  data,  a  eonsequenee  of  an  almost  eonstant  wave  height 
aeross  the  instrument  transeet.  The  equivalent  skill  value  implies  the  model  reproduees 
the  eross-shore  wave  transformation  for  Duok94  as  well  as  in  the  Delilah  ease. 
Correlation  values  from  all  beaehes  are  statistieally  signifieant  at  the  95  %  eonfidenee 
interval. 

Sensitivity  analyses  for  variations  in  7  were  performed.  It  was  found  that 
changes  in  y^^^  of  10%  result  in  a  mean  variation  of  only  3%  in  the  wave  model  skill  over 
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the  350  hour  span  while  ehanges  of  20%  showed  skill  ehanges  up  to  7%,  indieating  the 
model  is  relatively  insensitive  to  variations  in  y  values. 

Model  results  using  also  prediet  the  measured  wave  transformation 
well  (skill  >  0.95  )  with  eorrelation  eoeffieients  similar  to  those  found  in  the  model  runs 
using  y^p^.  One  exeeption  is  the  Torrey  Pines  ease  (skill  =0.80)  where  the  model  under 

prediets  wave  dissipation,  indieating  a  predieted  that  is  too  large.  The  ealeulated 
values  of  y^^gj  for  the  5  days  of  Torrey  Pines  analyzed  data  are  0.52-0.53,  whieh  are  50- 
100  %  larger  than  y^^,  values,  resulting  in  a  eorrelation  eoeffieient  of  0.74  and  a  skill  of 
0.76  (top  left  panel.  Figure  4.4).  The  larger  predieted  y  values  lead  to  less  dissipation 
aeross  the  surf  zone  (top  left  panel.  Figure  4.2)  resulting  in  a  bias  showing  the  model 
predieted  wave  height  being  larger  than  the  measured  wave  height.  Comparing  wave  and 
beaeh  eonditions  (Table  4.1)  between  the  planar  beaehes  provides  insight  into  the 
mismateh  between  y^^^  and  y^sis  allowing  for  an  extension  to  the  BS85  formulation.  The 

value  of  Sq  and  aeeordingly  and  are  similar  between  Torrey  Pines  and  Santa 

Barbara  resulting  in  equivalent  ineident  wave  energy,  initial  breaking  depths  and 
predieted  y^^gj  (Equation  18).  However,  the  beaeh  slope  at  Torrey  Pines  (1:50)  is  half 
that  of  Santa  Barbara  (1:25)  ereating  a  eross-shore  dissipation  distanee  twiee  that  of  Santa 
Barbara.  This  suggests  that  in  addition  to  wave  steepness  (5^) ,  beaeh  slope  (tany^) 
beeomes  important  in  the  parameterization  of  y  for  shallow  sloping  beaehes  with  long 
period  swell.  This  is  eonsistent  with  results  from  Durand  and  Allsop  (1997)  who 
examined  waves  in  a  flume  over  varying  bed  slopes  (1:10-1:30).  They  eoneluded  that 
wave  breaking  and  energy  dissipation  is  eompressed  into  shorter  eross-shore  distanees 
over  steeper  slopes  for  ineident  waves  with  similar  wave  steepness  values. 

For  wave  steepness  values  less  than  0.002,  y^^gj  values  no  longer  show 
dependenee  on  (Figure  4.1).  Some  data  from  Torrey  Pines,  Santa  Barbara  and  Duek94 
fall  into  this  range  and  are  separately  plotted  in  Figure  4.5.  At  these  small  values  of 
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there  appears  to  be  some  dependenee  on  the  Iribarren  number  (Battjes,  1974),  ^  = 

whieh  ineludes  the  beaeh  slope,  tan  /3 .  A  best-fit  hyperbolie  tangent  line  has  been  fitted 
to  these  data  deseribed  by 

y  =  0.2  + 0.32  tanh(^);  <  0.002  (20) 

The  beaeh  slope  was  defined  using  a  line  eonneeting  the  depth  offshore  at 
h  =  47/^^, 0  (outside  the  surf  zone)  to  the  mean  water  level,  ineluding  tide  variation, 
interseetion  at  the  beaeh  faee.  On  barred  beaehes,  onee  the  wave  height  derived  water 
depth  beeomes  less  than  the  bar  depth,  the  foreshore  slope  viee  the  bar  slope  is  used, 
whieh  was  the  ease  for  the  Duek94  data,  resulting  in  large  Iribarren  number  values. 

Breaking  wave  types  have  been  shown  to  be  a  funetion  of  ^ ,  whieh  are 
delineated  in  Figure  4.6.  Correlation  eoeffieients  (Figure  4.6)  show  the  improvement  in 
using  Equation  20  for  determining  y  for  wave  steepness  values  less  than  0.002.  In  all 
eases,  the  wave  model  did  as  well  or  better  than  y  values  determined  from  the  BS85 
parameterization  Equation  18.  In  the  ease  of  Torrey  Pines,  there  was  a  17  %  improvement 
in  the  eorrelation  eoeffieient.  The  predietion  of  breaker  type  using  ^  with  the  Torrey 
Pines  data  in  the  spill-plunge  region  and  Santa  Barbara  data  in  the  plunge  range 
eorresponds  well  with  observations  from  literature  (Thornton  and  Guza,  1983  and 
Thornton  and  Guza,  1986). 

2.  Bed  Roughness  Length  (kj 

a.  Best-fit  Values 

The  objeetives  of  this  seetion  are  to  evaluate  the  bed  roughness  length 
parameter  and  its  empirieal  relationship  to  the  bed  shear  stress  eoeffieient  , 

investigate  methods  for  determining  a  prior  and  to  aseertain  the  flow  model 
sensitivity  to  k^ .  Eor  steady  eurrents,  depends  only  on  bottom  roughness  when  the 
bed  is  hydraulieally  rough,  and  several  empirieal  relationships  are  available  throughout 


tan  [5 
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the  literature.  DelftSD  determines  the  bed  shear  stress  eoeffieient,  Cj- ,  as  the  inverse 

square  of  the  White-Colebrook  formulation  (hereafter  WC),  Equation  7.  WC  is 
proportional  to  the  square  of  the  ratio  of  the  bed  roughness  length  (or  equivalent 

Nikuradse  roughness  of  the  bed)  to  the  total  water  depth,  [kjh^ ,  implying  the  Cj 

formulation  in  the  model  is  not  eonstant  in  the  eross-shore  spatially  and  temporally  as  a 
result  of  varying  bathymetry  and  tidal  ehanges.  Literature  suggests  that  ean  be  related 
to  sediment  grain  size  or  the  bed  roughness  length,  two  length  seales  that  ean  be  variable 
in  time  and  spaee.  Best  fit  values  for  ehosen  to  ealibrate  and  observations  for 

the  entire  300-hour  Delilah  data  set  are  k^  =  0.003  m  (Roller  model).  The  corresponding 
cross-shore  averaged  =  0.002  and  0.003 .  The  Santa  Barbara  optimal  value  for 
k^  =  0.009  m  for  both  roller  and  non-roller  model  runs  results  in  a  cross-shore  averaged 
Cf  =  0.004 . 

b.  Observations 

On  the  barred  beach  at  Duck,  North  Carolina  during  the  Superduck 
experiment  in  1985,  Whitford  and  Thornton  (1996)  used  a  mobile  instrument  sled  to 
analyze  various  terms  in  the  momentum  equation.  They  found  bed  shear  stress 
coefficients  for  offshore  of  the  bar  (0.004),  on  top  of  the  bar  (0.003)  and  in  the  trough 
(0.001).  Their  formulation  of  was  based  on  a  residual  of  the  local  alongshore 

momentum  balance.  No  direct  measurements  of  the  bed  roughness  were  made. 

Direct  estimates  of  the  bed  shear  stress  within  the  surf  zone  were 
calculated  by  Garcez -Faria  et  al.,  (1998)  during  the  Duck94  experiment  from  logarithmic 
velocity  profiles  measured  with  a  vertical  array  of  current  meters  mounted  on  a  moveable 
sled.  Bed  roughness  was  measured  using  a  sonar  altimeter  mounted  on  an  amphibious 
vehicle.  The  values  were  found  proportional  to  the  bottom  roughness  and  varied  by  an 

order  of  magnitude  across  the  surf  zone  (0.0006-0.012).  The  values  were  also 
calculated  using  a  modified  version  of  the  Manning-Strickler  bed  shear  stress  formulation 
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where  is  the  apparent  bed  roughness  length  scale  used  to  represent  bed  roughness  in 

the  presence  of  waves.  Garcez-Faria  et  ah,  (1998)  concluded  from  a  comparison  between 
their  empirical  relation  and  data  that  a  single  roughness  length,  equivalent  to  the 
measured  roughness,  could  be  used  to  characterize  combined  wave-current  flows  over  a 
moveable  bed  at  least  for  the  case  of  strong  alongshore  currents. 

Gallagher  et  ah,  (2001)  similarly  measured  roughness  values  in  the  trough 
and  bar  region  during  Duck94  using  sonar  altimetry  techniques,  and  found  averaged 
alongshore  rms  roughness  heights  between  0.02-0.06  m,  depending  on  the  cross-shore 
location.  They  found  large  spatial  and  temporal  variability  of  the  bed  roughness  inside  the 
surf  zone  (water  depths  <  2  m)  with  bed  roughness  decreasing  in  the  offshore  direction. 

Over  a  planar  beach  during  the  NSTS  experiment  at  Santa  Barbara, 
Thornton  and  Guza  (1986)  calculated  cross-shore  average  values  based  on  a  least 

square  fit  to  the  alongshore  current  model  solution.  The  mean  values  for  the 

alongshore  current  inside  the  breaker  line  were  0.009  using  a  linearized  bed  shear  stress 
formulation,  and  =  0.006  using  a  more  realistic  non-linear  bed  stress  formulation.  The 

wave-breaking  model  utilized  by  Thornton  and  Guza  to  calculate  the  current  field  did  not 
include  a  roller  dissipation  mechanism. 

Another  approach  in  determining  would  be  to  use  a  formulation  related 
to  the  sand  grain  size.  One  often  applied  estimate  is  k^  =  2.5d^Q  (Soulsby,  1997),  which  in 
this  case  gives  k^  =  0.00045  m  for  Delilah  and  0.006  m  for  Santa  Barbara. 
c.  Model-Data  Comparisons 

A  wide  range  of  k^  values  is  calculated  based  on  the  measured  sediment 
grain  size  and  bed  roughness  values.  Since  k^  is  the  only  tuning  parameter  used  in  the 
flow  portion  of  the  model,  it  is  important  to  understand  the  model  sensitivity  to  k^  and 
the  resultant  effect  on  the  cross-shore  and  alongshore  current  profile.  Values,  one  order  of 
magnitude  on  either  side  of  best  fit  mean  k^  (0.003  m),  are  chosen  for  the  barred  beach 
case  based  on  sediment  grain  size  {k^  =0.00045  m)  and  average  rms  bed  roughness 
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height  {k^  =0.04  m)  resulting  in  a  range  using  the  WC  formulation  of  0.001-0.004, 

eonsistent  with  literature.  For  Santa  Barbara  no  bed  roughness  measurements  were  made. 
Instead,  the  value  (0.006)  found  by  Thornton  and  Guza,  (1986)  using  a  non-linear  bed 

shear  stress  formulation  will  be  used  as  the  upper  bound  ( =  0.04  m).  The  lower  bound 
is  based  on  the  measured  mean  sediment  grain  size  {k^=  0.0006  m). 

The  value  of  k^  influenees  the  shape  of  the  eross-shore  and  alongshore 
distribution  of  the  alongshore  eurrent  using  the  WC  formulation  for  .  The  degree  of 

influenee  is  investigated  in  sensitivity  tests  by  eomparing  model  output  with  observations 
from  the  low  and  high  tides  of  the  Delilah  (Get.  10,15)  and  the  high  tide  of  the  Santa 
Barbara  (Feb.  5)  experiments  (Figures  4.7  and  4.8).  In  general  for  the  barred  beaeh  ease, 
deereasing  k^  an  order  of  magnitude  (0.00045)  based  on  the  sediment  grain  size,  results 

in  an  inerease  in  and  from  the  bar  region  to  the  beaeh  faee,  whieh  may  (low  tide. 
Get.  10)  or  may  not  (low  tide  Get.  15)  improve  the  fit.  The  magnitude  of  ,  for  fixed  k^ 
shows  little  variation  in  the  eross-shore.  Inereasing  k^  an  order  of  magnitude,  to  the 
average  measured  bed  roughness  length  (0.04),  deereases  and  v),  resulting  in 
underestimates  of  the  eurrents.  For  the  larger  k^  value,  Cy  is  inversely  related  (through 
depth,  h)  to  the  eross-shore  bathymetry  profile.  Gn  Get.  15,  inereasing  k^  at  low  tide 
would  provide  a  better  mateh  to  observations.  This  eorresponds  to  inereased  bottom 
roughness  qualitatively  seen  in  the  trough  from  bathymetry  measurements  (see  eross- 
shore  profile.  Figure  4.8)  after  the  Get.  13  storm.  Likewise,  deereasing  k^  at  low  tide  on 
Got.  10,  due  to  deoreased  roughness  in  the  trough  would  prove  a  better  fit  (see  eross- 
shore  profile.  Figure  4.7).  The  smaller  k^  value  (0.00045  m)  resulted  in  a  smoother  eross- 

shore  eurrent  profile  with  a  olearly  defined  maximum.  The  larger  k^  (0.04  m)  oreated  a 

eross-shore  eurrent  profile  with  two  maxima  residing  at  the  deeper  portions  of  the  profile 
(trough  and  offshore  the  bar),  with  strongest  signal  at  low  tide.  Water  elevation  ehanges 
from  the  high  to  the  low  tide  at  Delilah  inereased  Cy  on  the  bar  erest,  with  the  greatest 

ehange  (0.001)  found  for  the  largest  value  of  k^  (0.04  m).  Due  to  the  WC  log  rhythmie 
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dependency  of  on  ,  order  of  magnitude  changes  in  result  in  smaller,  non-linear 
changes  in  . 

For  the  planar  beach  case,  using  k^  derived  from  the  sediment  grain  size 
relation  (0.0006  m),  shifts  the  location  of  maximum  to  the  beach  face  and  doubles  the 
magnitude  (Figure  4.9).  There  is  a  gradual  increase  of  in  the  cross-shore  from 
offshore  to  the  beach  face,  with  an  average  value  of  0.002.  A  spike  in  at  the  beach 
face  causes  to  abruptly  reduce  to  zero.  Increasing  k^  an  order  of  magnitude  from  best 
fit,  and  subsequently  ,  does  not  change  the  shape  of  the  cross-shore  profile  and  only 
reduces  the  magnitude  of  the  velocity  0(0.1  m/5  ). 
d.  Summary 

The  values  for  k^  that  provided  a  best  fit  to  i/.  and  v/  observations  for  the 
barred  beach  data  set  is  k^  =  0.003  m  and  for  the  planar  beach  is  k^  =  0.009  m.  These 
values  did  not  correspond  to  either  the  measured  bed  roughness  length  or  the  sediment 
grain  size  for  the  beaches  examined.  The  corresponding  cross-shore  averaged 

resulting  from  the  WC  bed  shear  stress  formulation  was  0.002  for  the  barred  beach  case. 
Values  of  are  similar  to  the  mean  results  found  by  Whitford  and  Thornton  (1993)  for 

the  same  barred  beach  when  using  estimates  of  various  terms  in  the  momentum  equations 
to  determine  the  bed  shear  stress  coefficient,  but  were  less  than  those  found  by  Garcez- 
Faria  et  al.  (1998)  using  a  variant  of  the  Manning-Strickler  Cy  formulation.  Planar  beach 

best- fit  values  for  k^  =  0.009  m  resulted  in  a  cross-shore  average  Cy  =  0.004  .  The  value 
for  Cy  is  slightly  less  than  that  found  by  Thornton  and  Guza  (1986)  by  fitting  an 
alongshore  current  model  with  data  using  a  non-linear  bed  stress  formulation. 

The  assumption  that  the  cross-shore  distribution  of  the  bottom  roughness 
length  scale  ( )  remains  unchanged  during  the  span  of  the  Delilah  experiment,  during 
which  significant  changes  in  wave  forcing  and  bathymetry  occurred,  is  not  realistic. 
Nevertheless,  a  single  k^  based  on  fitting  all  the  data  gives  reasonable  results  Oct.  7-19, 
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which  suggests  that  the  eross-shore  variation  of  is  not  overly  sensitive  to  ehanges  in 
and  is  mostly  eontrolled  by  bathymetric  changes  and  tidal  variation.  The  tests  indicate 
that  model  predietions  vary  with  variations  in  but  are  not  overly  sensitive  to  order 
magnitude  variations  in  .  A  bed  shear  stress  formulation  in  which  k^  represented  a 
measurable  quantity  would  be  more  appropriate. 


C.  1-DH  VERSUS  2-DH  MODELING 

1.  Introduction 

The  importanee  of  using  the  roller  dissipation  meehanism  and  inclusion  of 
pressure  gradients  indueed  by  non-uniform  bathymetry  (1-D  versus  2-DH  modeling)  is 
examined  on  both  barred  and  planar  beaches.  The  conditions  using  measured  bathymetry 
(MB)  and  uniform  bathymetry  (UB)  with  roller  (R)  and  no  roller  (NR)  wave  dissipation 
meehanisms  are  eompared  with  measured  data.  In  the  ease  of  the  barred  beaeh,  where 
eontinuous  wave  and  eurrent  observations  were  made,  the  model  output  are  further 
evaluated  at  low  and  high  tide  to  examine  modulation  of  the  cross-shore  radiation  stress 
forcing  of  the  alongshore  current  as  a  result  of  the  tide.  The  uniform  bathymetry  ease  (1- 
D)  is  produced  using  the  cross-shore  instrument  transect  bathymetry  profile  extended 
uniformly  in  the  alongshore. 

One  problem  with  comparing  computed  results  with  observations  is  determining 
the  extent  to  which  remote  bathymetry  can  influence  the  flow  field  in  the  model  domain. 
Remote  bathymetry  is  defined  as  strong  variations  in  the  alongshore  bathymetry  upstream 
of  the  measurement  mini-grid  area  and  therefore  not  ineluded  in  the  model  domain.  It  is 
believed  that  these  features  can  have  a  strong  effeet  on  the  hydrodynamies.  During  the 
Delilah  experiment,  time-lapse  video  exposures  were  used  to  reveal  shallow  regions  of 
bathymetry  (migrating  sand  bars)  not  eaptured  by  the  daily  bathymetry  measurements  in 
the  mini-grid  (delineated  by  lines  in  Figure  4.10  a.-m.).  Remote  bathymetry  structures 
can  be  seen  in  the  video  images  as  lightly  colored  areas  representing  wave  breaking 
regions  owing  to  deereased  depth.  On  Oet.  7-8,  strong  nearshore  sinusoidal  variation  in 
the  alongshore  was  present  during  small  waves  and  weak  currents  (Figure  4.10  a.-b.). 
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Waves  and  currents  increased  on  Oct.  9-13  (Figure  4.10  c.-g.)  from  the  south  due  to  a 
storm  front  and  distant  Hurricane  Lili  resulting  in  the  linear  nature  of  the  bar  and  a  bulge 
of  sand  appearing  south  of  the  mini-grid  lateral  boundary  on  Oct.  12  (Figure  4.10  f).  On 
Oct.  13  (Figure  4.10  1.),  the  bulge  of  sand  entered  the  mini-grid  region,  but  was  not 
measured  as  the  CRAB  (bathymetry  collection  vehicle)  was  unable  to  operate  in  the  high 
waves.  During  the  first  day  the  sand  bulge  was  measured  (Oct.l4)  (Figure  4.10  h.),  waves 
were  mild,  a  condition  which  continued  through  Oct.  19.  After  Oct.  15  incident  wave 
direction  was  the  result  of  a  storm  to  the  north  and  the  departing  hurricane  to  the  south, 
resulting  in  mean  currents  less  than  0.30 m/5  on  Oct  16  (Figure  4.10  j.)  and  a  bi-modal 
wave  field  Oct.  17-19  (Figure  4.10  k.-m.).  With  no  knowledge  of  the  extent  of  upstream 
remote  bathymetry  (Oct.  7-13),  the  weak  currents  (v),v/  <  0.30m/5)  on  Oct.  16  and  the 

bi-modal  nature  of  the  incident  wave  field  (Oct.  17-19),  Oct.  15  is  chosen  for  modeling 
and  comparison  at  the  cross-shore  and  alongshore  transects  with  computed  output. 

Santa  Barbara  wave  and  bathymetry  conditions  changed  little  during  Feb.  4-6.  No 
time-lapse  video  was  taken  during  the  experiment  and  no  knowledge  of  upstream  remote 
bathymetry  exists.  Feb.  4  is  chosen  for  comparison. 

The  calibrated  wave  model  results  from  Delilah  are  described  first  (Section  C2a), 
followed  by  the  alongshore  current  velocity  distributions,  in  both  the  cross-shore  ( 17  ) 

(Section  C2b)  and  alongshore  (i/,)  (Section  C2c)  directions,  based  on  forcing  obtained 

from  the  wave  model.  The  same  comparison  format  is  utilized  for  the  Santa  Barbara  case 
(Section  C3)  minus  the  alongshore  ( i/, )  analysis. 

2,  Barred  Beach 

a.  Introduction 

Model  output  using  both  uniform  and  measured  bathymetry  are  compared 
with  measured  data  from  Delilah  at  high  (0400)  and  low  (1100)  tide  on  Oct.  15  with  the 
objective  of  examining  the  cross-shore  wave  height  transformation  as  well  as  the  cross¬ 
shore  and  alongshore  variations  of  the  alongshore  current.  Model  assessment  over  a  wide 
range  of  conditions  is  also  examined  using  the  statistics  of  model  skill  (5),  linear 
correlation  (r)  and  absolute  rms  difference  {s),  between  observed  data  and  computed 
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output  at  each  functioning  cross-shore  and  alongshore  sensor  for  the  300  hour  span 
during  Delilah.  Observations  from  CM  10  are  not  included,  which  was  subject  to 
intermittent  flooding  and  drying  at  low  water  (not  present  in  the  numerical  modeling). 
Instrument  transects  are  not  included  in  which  the  maximum  mean  (one-hour)  observed 
v)  or  vj,  <  0.2  nijs  (in  view  of  signal  to  noise  ratio). 

Oct.  15  offshore  conditions  are  characterized  by  moderate  swell  waves, 
decreasing  {H^.^=\.\  m  to  0.96  m)  over  the  half  tidal  cycle,  with  a  constant  1^  =  10.7  s 

incident  from  the  southeast  at  an  angle  of -12°.  Within  the  mini-grid,  the  bathymetry  is 
characterized  as  a  barred  beach  having  a  continuous  trough  with  the  presence  of 
alongshore  non-uniformities  and  a  bar  crest  located  at  200  m  offshore  (Figure  4.11).  The 
bulge  of  sand  that  progressed  into  the  mini-grid  due  to  persistent  northerly  alongshore 
currents  is  located  at  150  m  offshore  between  alongshore  locations  700-800  m. 

b.  Cross-shore  Variation  of 

J  rms 

The  cross-shore  variation  of  at  high  and  low  tide  for  both  local  and 

roller  dissipation  mechanisms  are  given  in  Figure  4.12.  The  number  in  the  right  comer  of 
each  panel  indicates  the  model  skill  with  the  cross-shore  bottom  profile  given  as  a 
reference.  A  bar  is  located  near  cross-shore  location  200  m.  Water  elevation  due  to  the 
tide  influences  the  cross-shore  location  of  wave  dissipation.  At  high  tide  (Figure  4.12, 
lower  panels)  wave  dissipation  decreases  over  the  bar  allowing  smaller  waves  to  pass 
without  breaking,  resulting  in  dominant  wave  breaking  on  the  beach  face  and  locally 
smaller  mean  set-up  over  the  bar  and  within  the  trough.  At  low  tide  (Figure  4.12,  top 
panels),  the  effect  of  the  bar  becomes  more  pronounced;  waves  break  offshore  of  the  bar, 
which  results  in  larger  mean  set-up  over  the  bar  and  in  the  trough.  Overall,  the 
computed compares  well  with  the  measurements  (model  skill  >  0.90)  with  only  a 
slight  variation  in  model  skill  between  the  local  and  roller  dissipation.  Agreement 
between  roller  and  no  roller  skill  values  indicates  the  roller  does  little  to  change  the  rms 
wave  height  stmcture,  but  instead  functions  to  delay  the  transfer  of  momentum  in  the 
cross-shore.  The  trend  is  similar  for  the  uniform  bathymetry  case  and  the  plot  is  omitted. 
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c.  Cross-shore  Variation  of 

The  cross-shore  transects,  along  with  calculated  2-DH  alongshore 

pressure  gradients  using  both  roller  and  no  roller  dynamics  are  discussed  (Figures  4.13 
and  4.14).  Pressure  gradients  do  not  form  in  the  case  of  uniform  alongshore  bathymetry, 
and  therefore,  are  not  included  in  the  plots.  At  both  high  (0400)  and  low  tide  (1100)  the 
observed  velocity  maximum  of  the  alongshore  current  (*)  occurs  in  the  trough  (Figures 
4.15  and  4.16).  At  high  tide,  observed  velocities  in  the  trough  are  slightly  larger  as  a 
result  of  larger  incident  offshore  wave  height  at  this  time. 

For  the  case  of  alongshore  uniform  bathymetry  utilizing  both  a  roller  and 
local  wave  dissipation  (no  roller)  (upper  panel.  Figures  4.15  and  4.16),  predominant 
breaking  is  predicted  at  the  bar  and  shore  with  the  computed  velocity  maxima  occurring 
locally  at  the  bar  and  shore.  The  poorest  model  skill  is  at  low  tide  for  uniform  bathymetry 
(top  panel.  Figure  4.15).  For  the  no  roller  case,  the  wave  energy  is  dissipated  locally  at 
breaking  resulting  in  changes  of  momentum  flux  with  overestimation  of  the  alongshore 
current  velocities  at  the  bar  and  shore  and  an  underestimation  of  the  alongshore  current 
velocities  in  the  trough.  At  high  tide  (Figure  4.16),  fewer  waves  break  on  the  bar,  and  the 
velocity  maximum  for  the  local  dissipation  model  shifts  from  the  bar  to  the  shore.  The 
roller  version  shows  a  shoreward  shift  of  the  velocity  peak  with  an  increase  in  the  trough 
velocity.  With  no  pressure  gradients  contributing  to  the  alongshore  forcing,  the  velocity 
increase  in  the  trough  for  the  roller  case  is  attributed  to  the  spatial  lag  in  wave  dissipation, 
resulting  in  a  transfer  of  momentum  to  the  trough. 

Examination  of  the  local  wave  dissipation  (no  roller)  results  over 
measured  and  uniform  bathymetry  indicates  that  pressure  gradients  can  contribute  to 
alongshore  currents.  At  low  tide,  the  predictions  over  measured  bathymetry  still  show  a 
velocity  peak  at  the  bar  for  the  no  roller  case,  though  less  than  that  found  for  uniform 
bathymetry.  Current  speeds  in  the  trough  are  larger  than  those  found  for  a  uniform 
bathymetry  and  no  model  velocity  peak  is  found  at  the  shore.  Model  skill  is  doubled  by 
the  inclusion  of  measured  bathymetry  during  low  tide  and  shows  a  50  percent 
improvement  at  high  tide.  The  velocity  differences  at  the  bar  and  in  the  trough  are  the 
result  of  a  negative  pressure  gradient  acting  along  the  bar  and  a  local  positive  gradient  in 
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the  trough  acting  at  the  cross-shore  array  (Figures  4.13  and  4.14).  The  difference  at  the 
shore  is  due  to  a  weak  negative  pressure  gradient. 

In  comparing  measured  bathymetry  case  for  roller  and  no  roller  models 
(bottom  panel,  Figures  4.15  and  4.16),  the  primary  interest  is  the  computed  location  of 
the  current  maximum.  Including  the  roller  case  results  in  an  even  stronger  flow  in  the 
trough  than  the  measured  bathymetry  with  no  roller,  especially  during  low  tide.  This  can 
be  attributed  to  first,  the  inclusion  of  the  roller  allowing  for  a  spatial  lag  in  the  dissipation 
of  energy  resulting  in  greater  velocities  in  the  trough,  and  second,  the  presence  of  a 
stronger  positive  pressure  gradient  acting  between  CM20  and  CM30  allowing  for  the 
shift  of  the  velocity  peak  from  the  bar  region  to  the  trough  (see  Figure  4.15).  The  cross¬ 
shore  distribution  of  the  roller  computed  alongshore  current  profile  matches  the  measured 
distribution  well  {s  =  0.71  low  tide,  5  =  0.81,  high  tide).  Over-prediction  by  the  roller 
model  in  the  bar  region  is  the  result  of  utilizing  a  spatially  constant  .  Large  ripples  that 
form  in  the  trough  after  the  Oct.  12-13  storm  (see  profile.  Figure  4.15)  suggests  the  use  of 
a  larger  local  value,  which  would  reduce  the  cross-shore  current  profde  and  improve 
skill  (as  demonstrated  in  Section  4C). 

The  predicted  velocity  distribution  is  improved  by  using  the  roller  in 
conjunction  with  measured  bathymetry.  Without  the  addition  of  pressure  gradients  as  a 
forcing  term,  the  roller  acting  over  uniform  bathymetry  over-predicts  velocities  at  the 
shore  during  high  and  low  tide  and  predicts  a  velocity  maximum  at  mid-bar  during  low 
tide. 

Measured  and  computed  velocities  for  the  various  cases  of  roller  and 

no  roller,  and  uniform  and  measured  bathymetry,  for  all  sensors  for  the  300-hours  at 
Delilah  are  plotted  in  Figure  4.17.  Values  for  skill  (s),  linear  correlation  (r)  and  absolute 
rms  difference  {s)  are  located  in  the  upper  left  corner  of  each  panel  with  a  perfect 
correlation  line  between  measured  and  observed  velocities  provided  as  a  reference.  The 
spread  of  data  is  largest  for  >0.7m/5  which  occurred  during  Oct.  9-13.  The  largest 

improvement  in  predicting  alongshore  velocities  is  by  including  a  roller  dissipation 
mechanism.  Without  the  roller,  the  model  consistently  underestimates  the  flow  field. 
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There  is  little  difference  in  the  skill  between  NR  cases  over  either  uniform  or  measured 
bathymetry,  indicating  the  roller  to  be  the  more  important  of  the  two  forcing  mechanisms 
for  high-energy  events.  Not  including  remote  bathymetry  may  be  one  reason  why  no 
improvement  in  model  skill  is  realized  for  the  MB  case  when  compared  with  UB  for  both 
R  and  NR.  Another  explanation  is  that  during  high-energy  events  with  strong  alongshore 
currents,  the  bathymetry  tends  to  a  linear  bar  and  1-D,  first  order,  balances  dominate. 

Measured  velocities  in  the  range  of  0.0 -0.5  m/5  encompass  the  periods 
Oct..  7-8,  14-15  and  18.  In  the  R  (MB)  case,  a  group  of  outlying  data  (Oct.  8)  is  the  result 
of  strong  rhythmic  nearshore  bathymetry  at  low  tide  creating  a  model-developed 
circulation.  The  gyres  do  not  form  in  the  R  (UB)  case,  as  the  bathymetry  is  straight  and 
parallel  (Figure  4.18). 

Good  agreement  is  found  between  observed  and  computed  velocities  on 
Oct.  14  and  15  (Table  4.8).  On  Oct.  14  and  15  the  bulge  of  sand  is  within  the  mini-grid 
domain  and  is  intruding  into  the  trough  region  (see  Figures  4.10  h,  i  and  4.1 1),  suggesting 
that  the  measured  bathymetry  may  play  a  role  in  the  alongshore  current  structure.  Skill 
(5),  linear  correlation  (r)  and  absolute  rms  difference  {s)  are  compared  between 
measured  and  modeled  results  using  the  cross-shore  array  without  and  with  CM  10  (only 
when  submerged)  in  conjunction  with  the  alongshore  array  to  determine  the  influence  of 
the  sand  bulge  on  the  alongshore  flow  in  the  mini-grid  domain  (Table  4.8).  The  R  (MB), 
R  (UB)  and  NR  (MB)  cases  not  including  CM  10,  all  have  similar  skill  and  correlation 
values  with  best  results  provided  by  R  (MB),  suggesting  both  the  roller  and  measured 
bathymetry  are  important  in  this  case.  Including  CM  10  at  the  beach  face  supports  the 
findings,  with  the  NR  (MB)  case  improving  in  skill  over  the  R  (UB)  case. 

Negative  velocities  represent  the  days  when  the  incident  waves  were  from 
the  north,  resulting  in  an  alongshore  current  to  the  south  (Oct.  16  and  19).  All  cases, 
except  the  R  (MB),  underestimate  the  transition. 
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d.  Alongshore  Variation  of  v 


The  alongshore  variation  of  the  alongshore  current  (Figures  4.19  and  4.20) 
indicates  an  alternating  deceleration  and  acceleration  of  the  flow  during  low  and  high  tide 
between  location  800  to  1000.  The  flow  decelerates  between  alongshore  location  800-850 
then  accelerates  between  850-925,  decelerating  again  from  alongshore  location  925-975 
and  finally  accelerates  975-985.  Wave  forcing  using  both  a  roller  and  local  dissipation 
mechanism  with  measured  (pressure  gradients)  and  uniform  bathymetry  are  utilized  to 
explain  the  observed  results. 

For  the  case  of  utilizing  uniform  alongshore  bathymetry,  the  importance  of 
the  roller  mechanism  is  demonstrated  by  a  3  fold  increase  in  skill  at  low  tide  and  2  fold 
greater  at  high  tide  over  the  no  roller  model  alongshore  currents.  The  increased  velocity 
found  in  the  roller  case  is  a  result  of  the  roller  delaying  momentum  transfer  until  the 
trough.  With  UB  (no  pressure  gradients)  no  alongshore  acceleration  is  evident. 

For  the  case  using  actual  measured  bathymetry,  the  local  wave  dissipation 
model  demonstrates  the  influence  of  the  alongshore  pressure  gradients  in  the  alongshore 
current  velocity.  In  both  the  high  and  low  tide  case,  model  computed  flow  over  measured 
bathymetry  is  able  to  replicate  the  measured  alternating  deceleration  and  acceleration  in 
the  alongshore  flow;  a  function  of  the  positive  pressure  gradient  upstream  of  CMS  5  and 
negative  pressure  gradient  near  CM  32-33.  This  flow  pattern  in  alongshore  current 
velocity  due  to  pressure  gradients  is  most  evident  when  including  a  roller  (skill  =  0.85 
(low  tide)  and  0.90  (high  tide)). 

Correlation  plots  between  measured  and  modeled  for  the  300  hours 

again  shows  the  importance  of  the  roller  dissipation  mechanism  in  the  proper  modeling  of 
the  alongshore  current  profile.  Values  of  skill  and  linear  correlation  are  one-third  larger 
for  R  compared  with  NR  cases,  with  the  absolute  error  indicating  the  NR  cases 
underestimating  the  alongshore  flow  by  as  much  as  0.5  (Figure  4.21).  The  spatial  lag 
in  momentum  transfer  due  to  the  roller  results  in  stronger  alongshore  flow  shifted  into  the 
trough  and  a  mean  mismatch  between  observed  and  predicted  of  0.2  m/5  .  The  skill 
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improvement  between  MB  and  UB  in  the  NR  eases  indieates  bathymetry  does  play  a  roll 
in  the  alongshore  flow  field  strueture. 

3,  Planar  Beach 

a.  Introduction 

Model  outputs  using  both  uniform  and  measured  bathymetry  is  seleeted 
for  eomparison  with  measured  data  from  Santa  Barbara  at  high  tide  (1100)  on  Feb.  4.  The 
representative  ease  is  used  to  assess  model  performanee  and  determine  the  importanee  of 
ineluding  a  roller  dissipation  meehanism  and  measured  bathymetry  on  a  near-planar 
beaeh.  Comparison  is  made  between  the  eross-shore  transformation  of  as  well  as  the 

eross-shore  variation  of  the  alongshore  eurrent,  using  model  skill  as  measure  of  eomparison. 
Skill  (5),  linear  eorrelation  (r)  and  absolute  rms  differenee  ( ^  )  between  observed  data  and 
eomputed  output  at  eaeh  of  the  12  eross-shore  sensors  for  13  hours  at  Santa  Barbara  is 
examined. 

The  offshore  eonditions  on  Feb.  4  are  eharaeterized  as  narrow-banded, 
moderate  swell  waves  =  0.8  m)  over  the  high  tide  eyele  with  a  long  period  (?),= 

14.2  s)  ineident  from  the  northwest  at  an  angle  of  -8°. The  bathymetry  within  the 
measurement  domain  ean  be  deseribed  as  a  near-planar  beaeh  with  a  slight  bend  in  the 
alongshore  depth  eontours,  resulting  in  a  steeper  profile  at  alongshore  loeation  300.  The 
bathymetry  beeomes  straight  and  parallel  offshore  in  3  m  depth  (Figure  4.22). 

b.  Cross-shore  Variation  of  and  v) 

^rms  model  predietion  as  a  funetion  of  eross-shore  distanee  is  shown  in 

Figure  4.23.  The  waves  gradually  shoal  and  then  deerease  in  height  as  the  wave  breaks. 
Waves  start  to  break  just  offshore  the  wave  height  maximum  (Figure  23).  The  model 
eomputed  for  all  oases  shows  good  agreement  with  the  measured  data  with  model 
skill  >  0.92. 

The  measured  velooity  maximum  of  the  alongshore  eurrent  ooours  where 
the  gradient  in  wave  transformation  is  the  greatest,  about  midway  between  the  shoreline 
and  the  loeation  of  the  maximum  (Figure  4.24).  The  velooity  appears  to  deoay  more 

slowly  in  the  seaward  direotion  than  in  the  shoreward  direotion.  In  general,  model  results 
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show  a  good  fit  with  measured  data.  There  is  an  overestimation  of  the  alongshore  eurrent 
at  eross-shore  loeation  90  in  the  NR  ease  over  both  uniform  and  measured  bathymetry, 
resulting  skill  values  are  0.83-0.84.  The  R  ease  has  a  better  fit  at  eross-shore  loeation  90, 
but  underestimates  the  observed  eurrent  maximum. 

Measured  and  eomputed  veloeities  at  eaeh  sensor  for  the  13-hours  of 
data  (Feb.  4-6)  at  Santa  Barbara  are  plotted  in  Figure  4.25.  Values  for  skill  (s),  linear 
eorrelation  (r)  and  absolute  rms  differenee  {s)  are  loeated  in  the  upper  left  hand  side  of 
eaeh  panel  with  a  perfeet  eorrelation  line  between  measured  and  observed  provided  as  a 
referenee. 

Largest  spread  of  data  oeeurs  for  measured  <  0.3  .  This  portion  of 

the  plot  represents  the  flow  assoeiated  with  initial  breaking  in  the  outer  surf  zone  and 
beyond.  At  the  higher  veloeities  (v^>0.4m/5)  it  appears  that  both  the  NR  eases 
represent  the  observations  better.  There  is  little  differenee  in  skill  between  the  R  (MB), 
NR  (MB)  and  NR  (UB)  eases,  though  best  skill  is  found  for  the  R  (MB)  ease. 

4,  Summary 

The  predieted  2-DH  eross-shore  and  alongshore  variation  of  the  mean 
alongshore  eurrent  over  both  a  barred  and  near-planar  beaeh  are  eompared  with  field 
observations  to  quantify  the  relative  importanee  of  eontributions  from  the  roller  wave 
dissipation  meehanism  and  the  alongshore  pressure  gradients  assoeiated  with  measured 
bathymetry.  Results  show  that  for  a  planar  beaeh  (Santa  Barbara)  and  high-energy  events 
on  a  barred  beaeh  (Delilah)  1-D  dynamies  perform  as  well  as  2-DH  dynamies.  One 
hypothesis  of  why  1-D  works  as  well  as  2-DH  on  the  barred  beaeh,  are  first,  the  bar 
trough  region  beeomes  linear  (alongshore  uniform)  during  times  of  strong  alongshore 
eurrents  resulting  from  large  ineident  waves  and  second,  the  results  are  similar  when  not 
including  the  remote  bathymetry  in  the  trough  region.  When  the  remote  bathymetry  was 
included  in  the  mini-grid  domain,  it  was  shown  that  the  2-DH  model  was  an 
improvement  over  1-D,  by  including  the  pressure  gradients  due  to  non-uniform 
alongshore  bathymetry,  in  properly  modeling  both  the  cross-shore  and  alongshore 
variation  of  the  alongshore  current.  Pressure  gradients  induced  by  variable  breaking  over 
measured  bathymetry  shifted  the  computed  current  maximum  into  the  trough  when  a 
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roller  dissipation  mechanism  was  used.  The  inclusion  of  the  roller  was  needed  to  properly 
model  the  magnitude  of  the  cross-shore  distribution  of  the  alongshore  current.  The  largest 
discrepancies  between  model  predictions  and  observations  of  cross-shore  and  alongshore 
distributions  of  the  alongshore  current  occur  when  a  local  wave  dissipation  mechanism  is 
used  over  alongshore  uniform  bathymetry  (1-D)  case. 


D,  MOMENTUM  ANALYSIS  (BARRED  BEACH) 

1.  Uniform  vs  Measured  Bathymetry 


The  contributions  to  the  various  terms  of  cross-shore  and  alongshore  momentum 
balances  (Equation  4)  contributing  to  the  mean  flow  over  uniform  and  non-uniform 
model  bathymetry  are  examined  at  low  tide  on  Oct.  15  at  each  of  3  cross-shore  transects 
(Figure  4.26).  A  representative  cross-shore  transect  (transect  1)  at  low  tide  (Figure  4.27) 
is  used  for  discussion.  First,  the  model  results  for  cross-shore  momentum  utilizing 
uniform  and  measured  bathymetry  are  described,  followed  by  the  alongshore  momentum. 
Cross-shore  integrated  values  of  the  momentum  terms  are  calculated  to  discuss  order-one 
balances.  Positive  values  for  terms  (see  Figure  4.27)  act  to  decelerate  the  flow  for  cross¬ 
shore  momentum  changes  and  accelerate  the  flow  for  alongshore  momentum  changes.  To 
ascertain  if  the  sum  of  the  cross-shore  and  alongshore  momentum  terms  balance  to  zero, 
the  changes  in  momentum  flux  terms  are  brought  to  the  right  hand  side  of  Equation  4  and 
become  negative  valued  terms.  In  turn,  to  maintain  cross-shore  and  alongshore 
orientation,  those  terms  must  be  multiplied  by  a  negative  sign  when  discussing  their 
contribution  to  the  flow  field. 


For  the  cross-shore  momentum  utilizing  uniform  bathymetry,  the  primary  balance 
along  all  transects  at  all  tide  levels  is  between  the  cross-shore  wave  forcing  (F^)  and 

dfj 

hydrostatic  pressure  force  due  to  the  cross-shore  gradient  of  the  mean  setup  ( — ).  The 

dx 


convective  mixing  term  ( v  —  ),  associated  with  the  alongshore  change  in  u  momentum 

dy 

becomes  important  when  comparing  with  measured  bathymetry  as  accelerations  and 
decelerations  in  the  alongshore  flow  occur  due  to  alongshore  non-uniformity  in  the 
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bathymetry  and  resultant  pressure  gradients.  Along  transeet  1,  sinee  v  is  direeted 
alongshore  positive,  that  portion  of  the  eross-shore  transeet  in  whieh  the  eonveetive 
mixing  term  is  negative  (top  right  panel,  Figure  4.27)  implies  offshore  direeted  u 
momentum  is  inereasing  in  the  alongshore  eompensating  for  the  alongshore  deerease  in 
veloeity  (Figure  4.26). 

Cross-shore  integrated  values  of  the  eross-shore  momentum  show  that  for  both 
the  uniform  and  measured  bathymetry  eases  wave  foreing  is  balaneed  by  the  eross-shore 
pressure  gradient  and  is  several  orders  of  magnitude  larger  than  the  other  terms  (Table 
4.2  and  4.3). 

The  model  predieted  alongshore  momentum  balanee  for  uniform  bathymetry  is 
primarily  between  the  alongshore  direeted  wave  foreing  ( F), )  and  the  alongshore  bottom 

frietion  (r^  )  (bottom  left  panel.  Figure  4.27),  whieh  has  been  largely  demonstrated  using 

a  1-D  formulation  both  analytieally  (Bowen  (1969),  Longuet-Higgins  (1970),  Thornton 
(1970)  and  others)  and  quantitatively  eompared  with  field  data  (Thornton  and  Guza, 
1986,  Fedderson  et  ah,  1998).  At  the  eross-shore  bar  region  (210  m)  where  wave  foreing 
is  strongest,  eontributions  by  the  turbulent  mixing  term  {m,)  aet  in  eonjunetion  with  the 
bottom  frietion  to  provide  a  balanee.  In  the  eross-shore  trough  region  (150  m)  as  waves 
reform,  A,,  beeomes  small  and  the  turbulent  mixing  eompletes  the  balanee  with  the 

bottom  frietion.  Examination  of  the  alongshore  veloeity  pattern  (Figure  4.19)  reveals  no 
aeeeleration  or  deeeleration  in  the  flow. 

Cross-shore  integrated  values  of  the  alongshore  momentum  for  uniform 
bathymetry  are  listed  in  Table  4.4  and  support  findings  of  the  eross-shore  distribution 
shown  in  Figure  4.27.  The  representative  ease,  (low  tide  transeet  1)  shows  a  strong 
eorrelation  between  and  .  Although  makes  signifieant  eontributions  to  the  flow, 

eross-shore  sign  ehanges  in  result  in  small  integrated  values.  Pressure  gradient  values 

are  near  zero,  eonfirming  that  uniform  bathymetry  does  not  generate  alongshore 
variations  in  the  water  level.  The  mean  momentum  flux  terms  provide  little  eontribution 
to  the  flow  field. 
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All  the  alongshore  momentum  terms  eontribute  for  the  measured  bathymetry  ease 
(right  panels,  Figure  4.27).  In  the  eross-shore,  wave  foreing  dominates  at  the  bar  and 
shore,  becoming  near  zero  in  the  trough  as  waves  reform.  The  alongshore  pressure 
gradient  changes  sign  in  the  cross-shore;  first  it  opposes  wave  forcing  in  the  bar  region 
and  decelerates  flow,  then  becomes  positive  shoreward  in  the  trough  and  acts  to 
accelerate  the  alongshore  flow.  The  inertia  term,  which  appears  to  act  opposite  the 
pressure  gradient,  is  an  indication  of  alongshore  acceleration  (negative  values)  and 
deceleration  (positive  values).  In  this  case,  offshore  of  the  bar,  inertia  indicates  a  slight 
acceleration  in  the  flow  due  to  bar  induced  wave  breaking.  At  the  bar,  where  the  pressure 
gradient  acts  opposite  the  wave  forcing,  the  positive  inertia  slightly  lags  the  pressure 
gradient,  but  still  indicates  the  deceleration  in  the  flow.  In  the  trough,  large  negative 
inertia  reflects  the  current  maxima  where  a  strong  pressure  gradient  acts  in  the  direction 
of  wave  forcing.  Turbulent  mixing,  m, ,  transfers  maxima  in  v  momentum  in  the  trough 
region  to  areas  of  slower  v  both  in  the  onshore  and  offshore  directions.  Lateral  mixing  is 
most  important  near  the  current  jets  where  it  reduces  and  broadens  vj. .  The 

convective  mixing  term  results  in  offshore  advection  by  the  undertow  velocity  of  the 
slower  V  momentum  at  the  shore  in  the  direction  of  the  trough  and  faster  v  momentum 
in  the  trough  further  offshore  toward  the  bar.  Both  mixing  terms  act  to  diffuse  the 
alongshore  flow  resulting  in  a  more  cross-shore  uniform  current  structure  downstream  as 
the  trough  becomes  more  alongshore  uniform. 

Cross-shore  integration  of  each  of  the  alongshore  momentum  terms  (Table  4.5) 
again  shows  the  high  correlation  that  exists  between  and  with  the  magnitude  of  the 

integrated  values  varying  less  than  10  %  from  the  uniform  bathymetry  case.  There 

dfi 

appears  to  be  a  second  balance  between  the  gradient  in  the  alongshore  set-up  ( — -)  and 

dy 

dv 

the  alongshore  inertia  (v  — )  where  cross-shore  integrated  values  between  the  two  are 

dy 

equal  in  magnitude  but  opposite  in  sign.  In  this  case  the  magnitudes  of  the  pressure 
gradient  and  inertia  are  approximately  75  %  that  of  the  wave  forcing  and  bottom  friction. 
This  suggests  that  the  non-uniform  bathymetry  may  contribute  to  the  flow  field  via  set-up 
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induced  pressure  gradients,  both  in  veloeity  and  alongshore  current  profde.  Turbulent 
mixing  as  well  as  eonveetive  mixing  appear  small  in  magnitude,  again  the  result  of 
integrating  positive  and  negative  eontributions,  but  make  an  important  contribution  in 
spreading  veloeity  peaks  to  a  more  uniform  cross-shore  profile. 

2.  Alongshore  Variation  of  the  Total  Momentum 

The  alongshore  variation  of  the  incident  wave  field  and  resultant  flow  field  are 
examined  over  measured  bathymetry  to  examine  the  differenee  in  contribution  made  by 
each  of  the  terms  in  the  alongshore  momentum  equation.  At  low  tide,  wave  height  decay 
increases  over  the  shallower  water  depths  at  both  the  bar  (cross-shore  loeation  210  m,  see 
Figure  4.26)  and  sand  bulge  (eross-shore  loeation  160  m  between  alongshore  loeations 
700-800)  resulting  in  loeally  higher  set-up  of  the  mean  water  level.  Smaller  variations  in 
the  alongshore  bathymetry  (holes  in  the  bar  between  1025-1075  and  950-990)  result  in 
loeal  negative  gradients  in  the  mean  alongshore  pressure  gradients  (note  offshore  turning 
of  flow  at  hole  locations  in  the  bar).  Slinn  et  ah,  (2000)  found,  when  using  sinuous 
approximated  bathymetry,  that  descending  bathymetry  in  the  alongshore  results  in 
positive  pressure  gradients  with  loeal  flow  acceleration  and  aseending  bathymetry  in  the 
formation  of  negative  pressure  gradients  and  flow  deceleration,  ft  appears  that  when 
using  measured  bathymetry  rules  assoeiated  with  areas  where  positive  and  negative 
pressure  gradients  form  beeomes  more  complieated  as  eross-shore  variations  in  the 
bathymetry  become  important  as  well.  For  example,  for  bathymetry  descending  in  the 
trough  region,  alongshore  loeation  1000-1100  m,  Slinn  et  ah,  (2000)  would  suggest 
positive  pressure  gradients.  However,  holes  in  the  bar  allow  smaller  waves  to  pass 
without  breaking  resulting  in  smaller  set-up  over  the  bar  and  trough  region  and  the 
development  of  negative  pressure  gradients  (Figure  4.26). 

Large  non-uniform  features  in  the  alongshore  bathymetry,  such  as  the  sand  bulge, 
ean  influenee  the  veloeity  of  the  alongshore  flow  as  mueh  as  200  m  downstream,  whereas 
localized  pressure  gradients  are  the  result  of  immediate  response  to  small  changes  in  the 
bathymetry.  At  the  sand  bulge,  a  double  hump  cross-shore  eurrent  strueture  exists  (Figure 
4.26,  bottom  panel)  as  a  result  of  the  increased  set-up  and  subsequent  large  alongshore 
pressure  gradient  (Figure  4.26,  middle  panel).  Further  downstream,  at  loeation  985  where 
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the  cross-shore  instrument  transect  resides  (Figure  4.26),  the  influence  of  the  sand  bulge 
ceases  and  two  gaps  in  the  bar  (alongshore  locations  950-990  m  and  1025-1075  m) 
0(0.2)  m  deep  modify  the  pressure  gradient  pattern.  This  pressure  gradient  results  in  a 
single  maximum  of  the  alongshore  current  in  the  trough,  which  is  verified  in  the 
measured  data  (Figure  4.15). 

The  terms  from  the  time-averaged  x  and  y  -momentum  balances  (Equation  4) 
using  measured  bathymetry  are  plotted  at  alongshore  transects  1 ,  2  and  3  over  low  tide 
(Figure  4.28).  The  cross-shore  momentum  terms  reflect  the  large  signature  of  the  pressure 
gradient  balancing  the  cross-shore  wave  forcing.  The  flow  decelerates  in  the  alongshore 
trough  and  bar  region  due  to  negative  pressure  gradients  and  divergent  alongshore  flow 
as  a  result  of  the  bar  holes  (transect  2,3).  The  convective  mixing  term  has  a  positive  value 
and  indicates  that  offshore  directed  u  momentum  is  decreasing  in  the  alongshore  as  the 
flow  decelerates. 

For  the  alongshore  momentum,  the  relative  strengths  and  cross-shore  structure  of 
each  term  varies  with  alongshore  position.  The  alongshore  wave  forcing  and  pressure 
gradient  terms  are  clearly  not  uniform  in  the  alongshore,  reflecting  the  variability  of  the 
bathymetry.  Transects  1  and  2  show  the  alongshore  pressure  gradient  act  opposite  the 
wave  forcing  at  the  bar  and  in  conjunction  with  the  wave  forcing  from  mid-trough 
shoreward.  Strongest  flow  resides  in  the  trough  at  transects  1  and  2  with  flow 
accelerating  at  transect  1  while  still  under  the  influence  of  the  sand  bulge.  A  large 
negative  pressure  gradient  upstream  of  Transect  2  decelerates  the  cross-shore  velocity 
profde.  At  transect  3,  weak  negative  pressure  gradients  act  opposite  the  wave  forcing 
across  the  entire  cross-shore  slowing  the  flow  in  the  trough  and  creating  a  more  uniform 
cross-shore  profile.  The  inertia  term  reflects  the  deceleration  and  cross-shore 
redistribution  of  flow  from  trough  maximum  to  cross-shore  uniform  between  transect  1  to 
3.  The  convective  mixing  term  has  its  largest  signature  shoreward  of  the  trough  at 
transect  1,  at  the  bar  and  trough  at  transect  2  and  no  signature  at  transect  3,  which 
coincides  with  maximum  v  and  largest  positive  alongshore  pressure  gradients.  At  these 
locations,  convective  mixing  works  to  redistribute  the  flow  to  a  more  cross-shore  uniform 
structure  by  reducing  the  v  velocity  peaks  and  broadening  areas  of  slower  moving  flow. 
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E,  MOMENTUM  ANALYSIS  (PLANAR  BEACH) 

1,  Model  vs.  Measured  Results 

The  cross-shore  variation  of  the  cross-shore  and  alongshore  momentum  terms  are 
discussed  for  transect  1  at  high  tide  (see  Figure  4.29).  Positive  values  for  the  cross-shore 
momentum  terms  and  negative  values  for  the  alongshore  momentum  terms  act  to 
accelerate  the  flow  (Figure  4.29). 

The  dominant  balance  for  the  cross-shore  momentum  using  measured  bathymetry 
on  a  planar  beach  is  between  the  cross-shore  wave  forcing  and  the  pressure  gradient.  All 
others  terms  are  negligible.  Cross-shore  integrated  values  (Table  4.6)  support  findings 
from  Figure  4.29. 

For  the  alongshore  momentum,  all  terms  contribute  to  the  force  balance.  The 
primary  balance  is  between  the  wave  forcing  and  the  friction.  Unlike  the  non-uniform 
bathymetry  case  from  Delilah,  the  pressure  gradient  here  is  small  in  comparison  to  the 
wave  forcing.  The  pressure  gradient  acts  in  conjunction  with  the  forcing  at  cross-shore 
location  110  m  and  is  otherwise  near  zero,  or  acting  to  decelerate  the  flow  in  conjunction 
with  bottom  friction.  The  inertia  term  shows  slight  acceleration  at  the  cross-shore 
location  where  wave  forcing  and  pressure  gradient  act  in  the  same  direction,  otherwise  it 
shows  deceleration.  Turbulent  mixing  is  small  and  the  convective  mixing  is  advecting 
higher  v  momentum  offshore. 

At  transect  1,  the  dominant  balance  is  between  wave  forcing  and  bottom  friction 
(Table  4.7).  The  pressure  gradient  is  0(20%)  of  wave  forcing  as  a  result  of  only  minor 
variations  in  the  alongshore  bathymetry. 

2.  Alongshore  Variation  of  the  Total  Momentum 

The  alongshore  variation  of  the  wave  transformation  shows  a  slight  increase  in 
wave  breaking  between  alongshore  location  300-350  m  resulting  from  a  more  rapid 
decrease  of  the  bathymetry  (Figure  4.30).  The  result  of  this  small  variation  in  alongshore 
wave  breaking  is  a  pressure  gradient  pattern  where  a  positive  pressure  gradient  exists  at 
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alongshore  location  300-350  m  and  a  negative  pressure  gradient  between  450-500  m.  The 
velocity  field  under  the  influence  of  the  weak  negative  pressure  gradient  results  in  a 
current  profile  with  the  velocity  maximum  residing  near  the  shore.  In  contrast,  the  area  of 
the  flow  field  under  the  influence  of  the  positive  pressure  gradient  shows  a  shift  of  the 
velocity  maximum  in  the  offshore  direction  with  the  flow  speed  tapering  off  in  the 
onshore  and  offshore  direction.  The  influence  of  the  pressure  gradient  and  resultant 
deceleration  in  the  nearshore  flow  from  transect  3  to  transect  1  is  shown  in  Figure  4.31. 
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V.  SUMMARY  AND  CONCLUSIONS 


The  capabilities  of  the  2-DH  numerical  nearshore  wave  and  current  model  known 
as  DelftSD  were  assessed  by  comparing  model  output  with  observations  from  barred  and 
planar  beaches.  The  model  solves  the  energy  balance  equation  for  the  evolution  of  the 
wave  directional  spectrum  parameterized  in  frequency  space.  The  current  model, 
initialized  by  the  wave  output,  solves  the  depth-averaged  x-  and  y-momentum  flux 
equations  to  determine  the  flow  field. 

The  model  was  assessed  for  sensitivity  to  the  two  free  parameters  y  for  wave 
breaking  and  for  bottom  shear  stress.  The  calibration  formula  developed  by  Battjes 
and  Stive  (1985)  to  determine  the  breaking  wave  dissipation  parameter  y  as  a  function  of 
the  deep  water  wave  steepness,  ,  in  the  Battjes  and  Janssen  (1978)  wave  transformation 
model  was  verified.  Good  comparisons  of  the  predicted  and  measured  values  were 
obtained  on  both  barred  and  planar  beaches  using  y  from  Battjes  and  Stive  (1985)  for 
>0.002  with  a  model  skill  equal  to  0.91.  For  <0.002,  y  was  found  independent  of 
Sq  .  Using  long  period  swell  (small  ^g)  data  acquired  at  Torrey  Pines  and  Santa  Barbara, 
California  and  Duck,  North  Carolina,  a  new  parameterization  for  y  was  introduced  based 
on  the  Iribarren  number,  which  includes  the  beach  slope.  Improved  H predictions 
were  obtained  using  the  y  formulation  based  on  the  Iribarren  number  for  <0.002  with 
a  model  skill  equal  to  0.90. 

A  bed  shear  stress  formulation  in  which  represented  a  measurable  quantity  was 
sought.  Values  for  k^  that  provided  a  computed  best  fit  to  vj  and  vj,  observations  for  the 
entire  data  sets  were  k^  =  0.003  m  for  the  barred  beach  and  k^  =  0.009  m  for  the  planar 

beach.  However,  these  values  do  not  correspond  to  either  the  measured  bed  roughness 
height,  0(0.04),  or  the  sediment  grain  size,  0(0.0002  m)  for  the  barred  beach  examined. 
The  corresponding  cross-shore  averaged  value  resulting  from  the  White-Colebrook 

bed  shear  stress  formulation  was  0.002  for  the  barred  beach  case.  Values  of  matched 
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the  mean  results  found  by  Whitford  and  Thornton  (1993)  for  a  barred  beaeh  when  using 
estimates  of  various  terms  in  the  momentum  equations  to  determine  the  bed  shear  stress 
eoeffieient,  but  were  less  than  those  found  by  Gareez-Faria  et  al.,  (1998)  using  a  variant 
of  the  Manning-Striekler  formulation.  The  planar  beaeh  best-fit  value  for  =  0.009 

m  resulted  in  a  eross-shore  average  =  0.004 .  The  value  for  was  slightly  less  than 
that  found  by  Thornton  and  Guza  (1986)  using  a  non-linear  bed  stress  formulation. 

The  assumption  that  the  bottom  roughness  height  {k^)  is  eonstant  in  the  eross- 

shore  and  remains  unehanged  during  the  span  of  the  Delilah  experiment  is  not  realistie  as 
signifieant  ehanges  in  wave  foreing  and  bathymetry  oeeurred.  Tests  indieate  that  model 
predietions  are  not  sensitive  to  order  of  magnitude  variations  in  k^ .  Using  the  best  fit  k^ 

for  all  runs  Get.  7-19,  suggests  that  the  eross-shore  variation  of  is  not  overly  sensitive 
to  ehanges  in  k^ ,  and  is  mostly  eontrolled  by  depth  ehanges  assoeiated  with  tidal 
variation. 

The  predieted  2-DH  eross-shore  and  alongshore  variations  of  the  mean  alongshore 
eurrent  over  both  barred  and  planar  beaehes  were  eompared  with  field  observations  to 
quantify  the  relative  importanee  of  eontributions  from  the  roller  wave  dissipation 
meehanism  and  the  alongshore  pressure  gradients  assoeiated  with  measured  bathymetry. 
Results  showed  that  on  a  planar  beaeh  (Santa  Barbara)  and  during  high-energy  events  on 
a  barred  beaeh  (Delilah),  1-D  dynamies  perform  as  well  as  2-DH  dynamies  in  modeling 
of  the  alongshore  eurrent.  On  the  barred  beaeh,  two  hypothesizes  are  provided  to  explain 
the  similarity  in  performanee.  First,  the  bar  trough  region  tends  to  be  uniform  during 
times  of  strong  alongshore  eurrents  ereated  by  large  ineident  waves,  and  seeond,  the 
results  are  similar  when  not  ineluding  the  remote  bathymetry  in  the  trough  region.  When 
these  remote  bathymetrie  variations  propagated  into  the  measurement  domain  (Oet.  14- 
15),  it  was  shown  that  a  2-DH  model  that  ineluded  the  pressure  gradients  due  to  non- 
uniform  alongshore  bathymetry  is  an  improvement  over  1-D  hydrodynamies  in  modeling 
both  the  eross-shore  and  alongshore  variations  of  the  alongshore  eurrent.  The  roller 
tended  to  shift  the  eomputed  eurrent  maximum  into  the  trough  owing  to  a  shoreward 
transfer  of  momentum.  The  inelusion  of  the  roller  was  needed  to  properly  model  the 
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magnitude  of  the  eross-shore  distribution  of  the  alongshore  eurrent.  The  largest 
diserepaneies  between  model  predietions  and  observations  of  eross-shore  and  alongshore 
distributions  of  the  alongshore  eurrent  oeeur  when  a  loeal  wave  dissipation  (no  roller) 
meehanism  is  used  over  alongshore  uniform  bathymetry  (1-D)  ease. 

The  eross-shore  and  alongshore  variations  of  the  wave  heights  and  resultant  flow 
field  were  examined  over  uniform  and  measured  bathymetry  to  examine  the  differenee  in 
eontribution  made  by  eaeh  of  the  terms  in  the  eross-shore  and  alongshore  momentum 
equation  using  the  roller  dissipation  meehanism.  As  to  be  expeeted,  the  primary  balance 
in  the  cross-shore  was  between  the  cross-shore  wave  forcing  and  hydrostatic  pressure 
force  due  to  the  cross-shore  gradient  of  the  mean  set-up.  These  balance  terms  were 
several  orders  of  magnitude  larger  than  the  other  momentum  terms.  For  uniform 
bathymetry,  the  momentum  balance  in  the  alongshore  was  found  to  be  between  the  wave 
forcing  and  the  bottom  shear  stress.  Using  measured  non-uniform  alongshore  bathymetry 
resulted  in  differential  wave  breaking  and  concomitant  alongshore  variability  in  set-up 
profdes,  thus  creating  alongshore  surface  slopes  (pressure  gradients).  These  pressure 
gradients  primarily  balanced  by  the  alongshore  inertia  term,  were  found  to  be  order  one 
forcing  mechanisms  on  barred  beaches,  but  only  30  percent  as  strong  on  planar  beaches. 
The  change  in  the  cross-shore  directed  momentum  in  the  alongshore,  convective  mixing, 
became  important  in  the  measured  bathymetry  case  as  a  result  of  accelerations  and 
decelerations  in  the  alongshore  current. 
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Figure  4.1.  Optimized  values  of  /  versus  deepwater  wave  steepness,  ,  for  data  from 

Battjes  &  Stive  (1985),  Delilah,  Torrey  Pines,  Santa  Barbara  and  Duek94.  Solid  eurve  is 
Equation  18.  Dashed-dot  line  indieates  extent  of  Battjes  &  Stive  data  sets.  Dashed  line 
indieates  values  equal  to  0.002. 
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Delilah 


Duck  94 


Figure  4.2.  Selected  days  of  Torrey  Pines  (Nov.  10),  Santa  Barbara  (Feb.  4),  Delilah  (Oct. 
10)  and  Duck94  (Sept.  16)  for  model  prediction  of  compared  with  observations  (*). 
The  solid  line  represents  model  prediction  using  an  optimized  y  and  the  dashed  line 
represents  model  output  using  the  parameterized  y ,  Equation  18.  The  bottom  profde  is 
given  as  a  reference.  Skill  is  given  in  the  legend. 
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Torrey  Pines  (y^pj)  Santa  Barbara  (y^pj) 


Figure  4.3.  .  Measured  from  Torrey  Pines  (Nov.  10,14,18,20-21),  Santa  Barbara 
(Feb.  4-6),  Duck94  (Sept.  16)  and  Delilah  (Oct.  7-19)  plotted  against  the  model 
using  Values  of  skill  (s),  linear  correlation  (r)  and  absolute  error  are 

provided.  Units  for  absolute  error  are  m/s. 
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Torrey  Pines  (Tgsgs) 


Figure  4.4.  Measured  from  Torrey  Pines  (Nov.  10,14,18,20-21),  Santa  Barbara 
(Feb.  4-6),  Duck94  (Sept.  16)  and  Delilah  (Oct.  7-19)  plotted  against  the  model 
using  (Equation  18).  Values  of  skill  (5),  linear  correlation  (r)  and  absolute  error 
are  provided.  Units  for  absolute  error  are  m/s. 
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Figure  4.5.  Plot  of  optimized  y  versus  the  Iribarren  number  (^).  Solid  line  represents 
parameterization  given  by  Equation  20.  Dashed  line  represents  the  parameterization 
given  by  BS85  (Equation  18)  for  Sq  <  0.002  .  Dotted  lines  delineate  regions  of  spilling, 
plunging  and  eollapsing  type  waves. 
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Torrey  Pines  Santa  Barbara 


Figure  4.6.  Measured  from  Torrey  Pines  (Nov.  10,14,18,20  and  21),  Santa  Barbara 
(Feb.  4,5  and  6)  and  Duck94  (Sept.  16)  plotted  against  the  model  using  Equation 
20.  Delilah  comparison  not  included  (all  5o>  0.002).  Values  of  skill  (s),  linear 
correlation  (  r  )  and  absolute  error  ( ^  )  are  provided.  Units  for  absolute  error  are  m/s. 
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Figure  4.7.  Left  panels:  values  based  on  best  fit  (0.003  m),  bed  roughness  length 

(0.04  m)  and  sediment  grain  size  based  (0.00045  m)  during  high  and  low  tide  on  Oct. 
10  at  Delilah.  Middle  panels:  cross-shore  variation  of  v  for  3  k^  values.  Right  panels: 
alongshore  variation  of  v  for  3  k^  values. 
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High  Tide  (OCT  15) 
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Figure  4.8.  Left  panels:  values  based  on  best  fit  (0.003  m),  bed  roughness  length 

(0.04  m)  and  sediment  grain  size  based  (0.00045  m)  during  high  and  low  tide  on  Oct. 
15  at  Delilah.  Middle  panels:  cross-shore  variation  of  v  for  3  k^  values.  Right  panels: 
alongshore  variation  of  v  for  3  k^  values. 
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Santa  Barbara  (FEB  05) 


Santa  Barbara  (FEB  05) 


Figure  4.9.  Left  panel;  values  based  on  best  fit  (0.009  m),  bed  roughness  height 
(0.04  m)  and  sediment  grain  size  based  (0.0006  m)  during  high  tide  on  Feb.  5  at  Santa 
Barbara.  Right  Panel:  cross-shore  variation  of  v  using  3  k^  values,  for  both  roller  and  no 
roller  dissipation  mechanisms. 
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Figure  4.10  (a.-f.).  Time-lapse  video  exposures  taken  during  the  Delilah  experiment  (Oct. 
7-12).  Lightly  colored  areas  represent  wave  breaking  regions  owing  to  shallow  depth. 
The  date  of  the  picture  is  labeled  in  the  right  hand  corner  of  each  picture,  lateral 
boundaries  of  the  mini-grid  are  delineated  with  white  lines. 
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Figure  4.10  (g.-L).  Time-lapse  video  exposures  taken  during  the  Delilah  experiment  (Oct. 
13-18).  Lightly  colored  areas  represent  wave  breaking  regions  owing  to  shallow  depth. 
The  date  of  the  picture  is  labeled  in  the  right  hand  corner  of  each  picture,  lateral 
boundaries  of  the  mini-grid  are  delineated  with  white  lines. 
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Figure  4.10  (m.).  Time-lapse  video  exposures  taken  during  the  Delilah  experiment 
(Oct.  19).  Lightly  colored  areas  represent  wave  breaking  regions  owing  to  shallow  depth. 
The  date  of  the  picture  is  labeled  in  the  right  hand  comer  of  each  picture,  lateral 
boundaries  of  the  mini-grid  are  delineated  with  white  lines. 
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Figure  4.11.  Delilah  mini-grid  bathymetry  for  Oct.  15.  Selected  current  meters  labeled  for 
reference.  Dotted  lines  1-3  denotes  cross-shore  transects  examined  for  the  cross-shore 
and  alongshore  momentum  balance. 
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Figure  4.12.  Model  prediction  of  (solid  line)  compared  with  observations  (*)  for  the 

low  (1100)  and  high  (0400)  tide  of  Oct.  15  for  both  no  roller  and  roller  dissipation 
mechanisms.  Number  in  right  comer  represents  model  skill. 
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Figure  4.13.  Model  eomputed  alongshore  pressure  gradients  for  the  low  tide  at  Delilah  on 
Get.  15.  The  lighter  shades  represent  positive  pressure  gradients  that  aet  from  right  to  left. 
Dark  shades  are  negative  pressure  gradients  that  aet  left  to  right.  Cireles  represent 
loeations  of  sensors.  Symbols  R,  NR,  and  MB  stand  for  roller,  no  roller  and  measured 
bathymetry. 
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Alongshore  Pressure  Gradient  (High  Tide) 


Figure  4.14.  Model  eomputed  alongshore  pressure  gradients  for  the  high  tide  at  Delilah 
on  Get.  15.  The  lighter  shades  represent  positive  pressure  gradients  that  aet  from  right  to 
left.  Dark  shades  are  negative  pressure  gradients  that  aet  left  to  right.  Cireles  represent 
locations  of  sensors.  Symbols  R,  NR,  and  MB  stand  for  roller,  no  roller  and  measured 
bathymetry. 
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Case1:RvsNR(UB) 


Case  3:  R  vs  NR  (MB) 


Cross-shore  Distance  (m) 


Figure  4.15.  Measured  and  model  computed  cross-shore  variation  of  the  alongshore 
current  for  low  tide  (1100)  at  Delilah  on  Oct.  15.  The  solid  line  represents  roller 
dissipation  (R)  with  measured  bathymetry  (MB),  dashed  line  R  with  uniform  bathymetry 
(UB),  dash-dot  line  no  roller  (NR)  with  MB  and  dotted  line  NR  with  UB.  Numbers  in 
right  comer  represent  model  skill  values.  Measured  alongshore  velocity  are  indicted  by 
(*).  Bathymetry  along  instmment  transect  provided  as  a  reference 
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Casel:  RvsNR(UB) 


Case  3:  R  vs  NR  (MB) 


Cross-shore  Distance  (m) 


Figure  4.16.  Measured  and  model  computed  cross-shore  variation  of  the  alongshore 
current  for  the  high  tide  (0400)  at  Delilah  on  Oct.  15.  The  solid  line  represents  roller 
dissipation  (R)  with  measured  bathymetry  (MB),  dashed  line  R  with  uniform  bathymetry 
(UB),  dash-dot  line  no  roller  (NR)  with  MB  and  dotted  line  NR  with  UB.  Numbers  in 
right  comer  represent  model  skill  values.  Instmment  bathymetry  transect  provided  as  a 
reference 
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Figure  4.17.  Measured  cross-shore  variation  of  the  alongshore  velocity,  v^from  Delilah 
(Oct.  7-19)  plotted  versus  model  v^.  Letter  combinations  in  the  title  correspond  to  the 
following:  Roller  dissipation  (R),  No  Roller  (NR),  Measured  Bathymetry  (MB),  Uniform 
Bathymetry  (UB).  Numbers  in  left  corner  represent  the  model  skill  (  5  ),  linear  correlation 
( r  )  and  absolute  error  (s). 
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Figure  4.18.  Oct.  8  model  output  during  low  tide  for  both  R(MB)  (left  panels)  and  R(UB) 
(right  panels).  The  dots  in  the  upper  panels  represent  locations  of  the  cross-shore  and 
alongshore  instruments.  Bottom  panels  are  the  corresponding  cross-shore  transects  with 
measured  data  (*)  plotted  versus  model  output  (solid  line).  Model  skill  (5)  and  cross¬ 
shore  profde  provided  as  a  reference. 
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Case  1:  R  vs  NR  (UB) 


Case  3:  R  vs  NR  (MB) 


Figure  4.19.  Measured  and  model  computed  alongshore  variation  of  the  alongshore 
current  for  the  low  tide  (1100)  at  Delilah  on  Oct.  15.  The  solid  line  represents  roller 
dissipation  (R)  with  measured  bathymetry  (MB),  dashed  line  R  with  uniform  bathymetry 
(UB),  dash-dot  line  no  roller  (NR)  with  MB  and  dotted  line  NR  with  UB.  Numbers  in 
right  corner  represent  model  skill  values. 


64 


Case1:  RvsNR(UB) 


Case  3:  R  vs  NR  (MB) 


Alongshore  Distance  (m) 


Figure  4.20.  Measured  and  model  computed  alongshore  variation  of  the  alongshore 
current  for  the  high  tide  (0400)  at  Delilah  on  Oct.  15.  The  solid  line  represents  roller 
dissipation  (R)  with  measured  bathymetry  (MB),  dashed  line  R  with  uniform  bathymetry 
(UB),  dash-dot  line  no  roller  (NR)  with  MB  and  dotted  line  NR  with  UB.  Numbers  in 
right  comer  represent  model  skill  values. 
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Figure  4.21.  Measured  alongshore  variation  of  the  alongshore  velocity,  v,,  from  Delilah 

(Oct.  7-19)  plotted  versus  model  v^. .  Letter  combinations  in  the  title  correspond  to  the 

following;  Roller  dissipation  (R),  No  Roller  (NR),  Measured  Bathymetry  (MB),  Uniform 
Bathymetry  (UB).  Numbers  in  left  comer  represent  model  skill  (5  ),  linear  correlation  (r  ) 
and  absolute  error  (  ^  ). 
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Figure  4.22.  Santa  Barbara  bathymetry  for  Feb.  4  .  Seleeted  eurrent  meters  labeled  for 
referenee.  Dotted  lines  denote  eross-shore  transeets  for  the  eross-shore  and  alongshore 
momentum  balanees. 
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Figure  4.23.  Left  panel;  Model  prediction  of  (solid  line)  using  measured  bathymetry 
compared  with  observations  (*)  for  the  high  tide  of  Feb  4.  Right  panel;  The  cross-shore 
variation  of  the  alongshore  current.  The  number  in  the  lower  left  hand  comer  indicates 
the  model  skill  level.  The  beach  profile  is  shown  as  a  reference. 
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Case  1:  R  vs  NR  (UB)  Case  2:  NR  vs  NR  (MB/UB) 


Case  3:  R  vs  NR  (MB) 
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Figure  4.24.  Measured  and  model  computed  cross-shore  variation  of  the  alongshore 
current  for  the  high  tide  (1100)  at  Santa  Barbara  on  Feb  4.  The  solid  line  represents  roller 
dissipation  (R)  with  measured  bathymetry  (MB),  dashed  line  R  with  uniform  bathymetry 
(UB),  dash-dot  line  no  roller  (NR)  with  MB  and  dotted  line  NR  with  UB.  Numbers  in 
right  comer  represent  model  skill  values.  Instmment  bathymetry  transect  provided  as  a 
reference. 
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Figure  4.25.  Measured  cross-shore  variation  of  the  alongshore  velocity,  it  from  Santa 
Barbara  (Feb.  4-6)  plotted  versus  model  .  Letter  combinations  in  the  title  correspond  to 
the  following:  Roller  dissipation  (R),  No  Roller  (NR),  Measured  Bathymetry  (MB), 
Uniform  Bathymetry  (UB).  Numbers  in  left  comer  represent  the  model  skill  (5),  linear 
correlation  ( r  )  and  absolute  error  (s). 
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Figure  4.26.  Top  panel:  Delilah  mini-grid  bathymetry  for  Get.  15  at  low  tide  overlaid 
with//^^^  eontours.  Middle  panel:  model  computed  pressure  gradients  at  low  tide.  The 
lighter  shades  represent  positive  pressure  gradients  that  act  in  the  direction  of  wave 
forcing  (right  to  left).  Dark  shades  are  negative  pressure  gradients  that  act  left  to  right. 
Bottom  panel:  Mean  (1  hour)  velocities  at  low  tide.  The  shaded  background  represents 
velocity  magnitude  in  m/s. 
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Uniform  Bathymetry  (Low  Tide) 


Measured  Bathymetry  (Low  Tide) 
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Measured  Bathymetry  (Low  Tide) 


Figure  4.27.  Top  panels;  model  predicted  cross-shore  variation  of  the  x-momentum  terms 
Equation  4  for  both  uniform  and  measured  bathymetry  for  the  low  tide  (1100)  of  Oct.  15 
transect  1  at  Delilah.  Bottom  panels:  cross-shore  variation  of  the  y-momentum  terms 
Equation  4  for  both  uniform  and  measured  bathymetry.  Model  depth  profile  provided  for 
reference. 
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Figure  4.28.  Top  panel:  model  computed  mean  alongshore  pressure  gradients  overlaid 
with  the  mean  velocity  computed  for  Delilah  on  Oct.  15  at  low  tide.  Middle  panel:  the 
cross-shore  variation  of  the  cross-shore  momentum.  Bottom  panel:  the  cross-shore 
variation  of  the  alongshore  momentum.  Transect  lines  numbered  1-3  are  provided  for 
reference. 
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Measured  Bathymetry  (High  Tide)  Measured  Bathymetry  (High  Tide) 


Figure  4.29.  Left  panel;  model  predicted  cross-shore  variation  of  the  x-momentum  terms 
Equation  4  for  measured  bathymetry  for  the  high  tide  (1100)  of  Feb.  4  transect  1  at  Santa 
Barbara.  Right  panel:  cross-shore  variation  of  the  y-momentum  terms  Equation  4  for 
measured  bathymetry.  Model  depth  profile  provided  for  reference. 
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Figure  4.30.  Top  panel:  Santa  Barbara  bathymetry  for  Feb.  4  at  high  tide  overlaid 
with//^^^  eontours.  Middle  panel:  model  computed  pressure  gradients  at  high  tide.  The 

lighter  shades  represent  positive  pressure  gradients  that  act  opposite  the  direction  of  wave 
forcing  (left  to  right).  Dark  shades  are  negative  pressure  gradients  that  act  right  to  left. 
Bottom  panel:  Mean  (1  hour)  velocities  at  high  tide.  The  shaded  background  represents 
velocity  magnitude  in  m/s. 
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Figure  4.31.  Top  panel;  model  computed  mean  alongshore  pressure  gradients  overlaid 
with  the  mean  velocity  computed  for  Santa  Barbara  on  Feb.  4  at  high  tide.  Middle  panel; 
the  cross-shore  variation  of  the  cross-shore  momentum.  Bottom  panel;  the  cross-shore 
variation  of  the  alongshore  momentum.  Transect  lines  numbered  1-3  are  provided  for 
reference. 
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Day 

■^0 

tan  P 

Torrey  Pines 

Nov.  10,14,18,20,21 

0.47-0.69 

0.06-0.08 

0.0011-0.002 

0.019-0.025 

Santa  Barbara 

Feb.  4-6 

0.29-0.65 

0.07-0.09 

0.001-0.002 

0.038-0.043 

Delilah 

Oct.  7-19 

0.33-1.76 

0.08-0.18 

0.003-0.02 

0.035  (off) 

Duck94 

Sept.  16 

0.28-0.33 

0.07 

0.0009-0.001 

0.066-0.105  (fore) 

Table  4.1.  Wave  and  beaeh  eonditions. 
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Transect 

F. 

< 

dn 

du 
u  — 
dx 

m, 

du 

V  — 

dy 

Residual 

High  Tide 

1 

-472.5 

1.3 

477.9 

0.3 

3.0 

-2.7 

7.2 

2 

-472.3 

1.2 

474.4 

0.2 

2.9 

-0.3 

6.1 

3 

-476.0 

1.3 

478.7 

3.5 

-1.6 

6.5 

Mid  Tide 

1 

-461.0 

2.2 

478.2 

1.1 

0.3 

-8.3 

12.5 

2 

-461.9 

1.7 

478.0 

0.8 

0.7 

-4.0 

15.3 

3 

-458.1 

1.7 

457.6 

2.4 

-0.4 

12.0 

15.1 

Low  Tide 

1 

-513.4 

2.9 

528.0 

-0.1 

0.2 

-2.3 

15.3 

2 

-512.1 

2.6 

526.0 

-0.1 

-0.1 

0.0 

16.3 

3 

-509.9 

2.5 

524.3 

-0.1 

-0.1 

-0.5 

16.3 

Table  4.2.  Cross-shore  integrated  cross-shore  momentum  values  using  uniform 
bathymetry  for  high,  mid  and  low  tides  at  each  of  the  alongshore  transects  for  Delilah  on 
Oct.  15.  Units  are  in  N/m. 
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Transect 

F. 

< 

dll 

du 
u  — 
dx 

m, 

du 

V  — 

dy 

Residual 

High  Tide 

1 

-479.7 

-7.0 

574.6 

15.8 

-4.6 

-83.0 

16.1 

2 

-621.9 

-10.2 

354.9 

2.4 

0.6 

288.1 

13.9 

3 

-578.3 

-2.2 

397.2 

3.4 

13.1 

180.3 

13.7 

Mid  Tide 

1 

-517.5 

-4.7 

592.6 

4.3 

-4.1 

-53.5 

17.1 

2 

-688.6 

-12.5 

407.5 

4.8 

-1.1 

309.5 

19.6 

3 

-601.8 

-0.9 

445.0 

2.2 

11.4 

159.7 

15.5 

Low  Tide 

1 

-551.6 

-2.4 

617.1 

3.4 

-4.0 

-44.0 

18.5 

2 

-708.2 

-11.2 

443.4 

2.0 

1.4 

293.2 

20.6 

3 

-612.9 

1.5 

447.9 

0.2 

9.4 

165.2 

11.3 

Table  4.3.  Cross-shore  integrated  cross-shore  momentum  values  using  measured 
bathymetry  for  high,  mid  and  low  tides  at  each  of  the  alongshore  transects  for  Delilah  on 
Oct.  15.  Units  are  InN/m. 
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Transect 

Fy 

dn 

dv 

V  — 

dy 

m, 

dv 
u  — 
dx 

Residual 

High  Tide 

1 

89.2 

-67.3 

-1.5 

-10.2 

-5.7 

-10.4 

-5.8 

2 

89.8 

-69.3 

0.3 

-10.2 

-7.0 

-9.7 

-6.2 

3 

89.9 

-70.7 

-2.9 

-8.1 

-5.1 

-9.3 

-6.4 

Mid  Tide 

1 

105.9 

-93.2 

-2.0 

-12.9 

14.8 

-7.8 

4.8 

2 

106.4 

-96.5 

2.4 

-16.1 

13.4 

-4.0 

5.5 

3 

105.6 

-98.8 

3.2 

-10.4 

8.5 

-1.8 

6.4 

Low  Tide 

1 

82.0 

-75.9 

0.6 

-10.1 

7.6 

-8.2 

-4.0 

2 

82.0 

-78.5 

1.3 

-9.2 

6.6 

-6.3 

-4.0 

3 

81.8 

-80.6 

1.8 

-7.3 

5.6 

-5.3 

-4.0 

Table  4.4.  Cross-shore  integrated  alongshore  momentum  values  using  uniform 
bathymetry  for  high,  mid  and  low  tides  at  each  of  the  alongshore  transects  for  Delilah  on 
Oct.  15.  Units  are  InN/m. 
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Transect 

Fy 

dn 

dv 

V  — 

dy 

m, 

dv 
u  — 
dx 

Residual 

High  Tide 

1 

75.7 

-81.8 

50.1 

-55.9 

1.2 

9.3 

-1.4 

2 

81.7 

-79.3 

123.7 

-124.2 

-61.9 

61.9 

-5.3 

3 

78.8 

-73.6 

-32.5 

20.6 

-6.6 

11.3 

-1.9 

Mid  Tide 

1 

101.4 

-100.8 

56.8 

-61.8 

2.5 

6.9 

5.0 

2 

117.1 

-99.7 

166.9 

-182.3 

-59.3 

60.1 

2.8 

3 

98.2 

-92.0 

-101.8 

85.4 

3.9 

11.8 

5.6 

Low  Tide 

1 

79.5 

-89.8 

54.2 

-42.6 

-1.4 

-3.1 

-3.2 

2 

91.4 

-85.8 

107.2 

-110.9 

-52.7 

46.3 

-4.5 

3 

68.0 

-74.4 

-78.3 

78.0 

1.4 

2.3 

-3.1 

Table  4.5.  Cross-shore  integrated  alongshore  momentum  values  using  measured 
bathymetry  for  high,  mid  and  low  tides  at  each  of  the  alongshore  transects  for  Delilah  on 
Oct.  15.  Units  are  InN/m. 


81 


Transect 

F. 

< 

dn 

du 
u  — 
dx 

m, 

du 

V  — 

dy 

Residual 

High  Tide 

1 

-1578.4 

-6.3 

1664.7 

11.2 

-15.4 

13.2 

89.1 

2 

-1690.3 

13.8 

1814.6 

0.7 

3.7 

-49.4 

93.1 

3 

-1657.3 

17.5 

1750.9 

0.0 

2.3 

-18.6 

94.8 

Table  4.6.  Cross-shore  integrated  eross-shore  momentum  values  over  measured 
bathymetry  at  high  tide  for  each  of  the  alongshore  transects  at  Santa  Barbara  on  Feb.  4. 
Units  are  in  N/m. 
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Transect 

u 

dll 

dv 

V  — 

dy 

dv 
u  — 
dx 

Residual 

High  Tide 

1 

-142.9 

118.4 

20.8 

30.0 

19.1 

-34.8 

10.5 

2 

-151.4 

141.3 

35.2 

-0.7 

-7.0 

-9.4 

8.0 

3 

-152.4 

140.4 

-61.1 

48.5 

8.3 

25.2 

8.9 

Table  4.7.  Cross-shore  integrated  alongshore  momentum  values  over  measured 
bathymetry  at  high  tide  for  eaeh  of  the  alongshore  transeets  at  Santa  Barbara  on  Feb.  4. 
Units  are  in  N/m. 
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Oct.  14-15 

R(MB) 

NR(MB) 

R(UB) 

NR(UB) 

v^(cm20-80) 

s 

0.67 

0.57 

0.64 

0.37 

r 

0.86 

0.76 

0.82 

0.45 

s 

0.12 

0.16 

0.13 

0.23 

(cm  10-80) 

s 

0.61 

0.58 

0.51 

0.36 

r 

0.80 

0.74 

0.67 

0.32 

s 

0.14 

0.15 

0.17 

0.22 

s 

0.75 

0.47 

0.69 

0.26 

r 

0.79 

0.26 

0.70 

0.38 

s 

0.10 

0.22 

0.13 

0.31 

Table  4.8.  Skill  (s),  correlation  (r)  and  absolute  rms  difference  (^)  values  for  and 

from  Delilah  (Oct.  14-15).  Symbols  R,  NR,  and  MB  stand  for  roller,  no  roller  and 
measured  bathymetry. 
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